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FREQUENCY  DEPENDENT  CHARACTERISTICS  IMPEDANCE. 
ATTENUATION.  AND  OHMIC  LOSS  OF  MICROSTRIP.  BASED  UPON 
NUMERICAL  CALCULATIONS  OF  EFFECTIVE  INDUCTANCE 

I.  INTRODUCTION 

\  general  matrix  method.  Miilahlc  for  efficient  numerical  solution  on  a  hig:h-*peed  com- 
[niter,  is  adapted  anil  applied  to  determining:  the  enrrent  distribution  and  trans\erse 
ma "net ie  field  of  uniform  mierostrip  transmission  lines.  Line  los*  and  freqiieney  are  in- 
eluded  in  the  solution  as  necessary  to  produee  accurate  results  for  the  current  distrihuliou. 
The  current  distrihuliou  is  then  used  to  provide  the  effective  characteristic  impedance,  in- 
duetance  per  unit  lengRli.  and  attenuation  of  the  line. 

The  characteristic  impedance  of  Microstrip  Transmission  Lines  has  been  of  interest  for 
25  years.  Young:1  provides  an  excellent  com|>ilation  of  selected  [tapers  including:  several 
early  papers  specifically  addressing:  Microstrip.  The  papers  hv  Colin2.  Wheeler.3  and 
Itryant  and  Weiss4  are  especially  fundamental  and  useful  for  common  engineerin': 
problems.  These  papers,  and  a  multitude  of  others  published  since,  depend  on  solving!  for 
the  static  capacitance  per  unit  len"th  of  the  selected  line  confi"tiratinn.  Most  practical 
geometries  do  not  lend  themselves  to  an  exact  analytic  solution  so  much  effort  has  been 
devoted  to  developing:  approximate  analytic  solutions.  With  the  advent  of  the  hig:h-speed 
computer  a  considerable  effort  has  been  devoted  to  developing:  numerical  techniques  for 
solving:  Laplace's  equation  to  yield  the  electric  field  configuration  and  capacitance  of 
useful  geometries.  Cui/an  and  Carver*  gtive  a  notable  example  of  this  approach. 

Once  the  capacitance  per  unit  lcng;th  has  been  determined  the  microwave  characteristic 
impedance  is  determined  hv  the  following:  relationship 

7.o  =  I  Mi  (It 


Leo  Young.  "Farallel  Coupled  Lines  and  Directional  Couplers.”  Artecli  House.  Inc..  1972. 

Y 

S.  B.  Cohn.  “Shielded  Coupled-strip  Transmission  Line.**  Il  l  I  Trans  on  Microwave  Theory  and  Techniques. 
Vol  MTT-3.  No.  5.  pp  29-38.  Oci  55. 

^  II.  A.  Wheeler.  'Transmission  line  Properties  of  Parallel  Strips  Separated  by  a  Dielectric  Sheet,”  II  I  I  Trans, 
on  Microwave  Theory  and  Techniques,  Vol.  MIT-13,  No.  2.  pp  I  72-185,  Mar  65. 

4 

1.  (».  Bryant  and  .1.  A.  Weiss.  “Parameters  of  Mierostrip  Transmission  lines  and  of  Coupled  Pairs  of  Micro¬ 
strip  Lines.”  Il  l  I  Trans,  on  Microwave  Theory  and  Techniques.  Vol.  MIT-16.  No.  12.  pp  1021*  1027.  Dec  68. 

**  O.  K  Cui/an  and  K.  V.  Carver,  “Characteristic  Impedance  of  Rectangular  Coaxial  Transmission  l  ines,”  II  I  I 
Trans,  on  Microwave  Theory  and  Techniques,  M  CI -12,  pp  489-495,  Sep  64. 
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" here  v  is  the  velocity  ot  propagation  in  the  lini'.  Vt  heeler1  anil  Brvant  ami  Weis:-1  ad¬ 
dress  tin1  practical  prnlilciu  of  determinin''  the  effective  velocity  in  a  line  where  the  fields 
are  partially  in  a  ilieleetrie  anil  partially  in  free  space.  The  results  shown  in  Figure  4  of 
Bryant  anil  Weiss1  will  he  iiseil  for  later  calculation  in  this  report  where  effective  velocity 
of  propagation  must  he  estimateil. 

The  approach  uscil  in  this  report,  while  nouanalvtic.  i-  general  anil  without  auv 
geometric  limitation  in  the  transverse  plane.  This  approach  is  uniipie  in  that  the  effects  of 
finite  comhicior  losses  anil  frequency  ilcpemlence  can  he  incluileil  in  the  analysis. 
Assumptions  maile  hv  other  authors  are  also  maile  here.  The  lines  to  he  analyzed  are 
assumeil  to  he  relatively  low  loss  lines  supporting  Ouasi-TFM  moiles.  (iapaeitance-haseil 
solutions  are  static  solutions  which  approach  exactness  only  for  lossless  lines.  'Pile 
i»»ihirlanre-ha»eil  solution  to  he  ilevelopeil  anil  applieil  here  can  allow  for  losses  hut  is 
quasi-static  anil  rctanlcil  potentials  have  not  heen  eonsiilereil. 

The  question  is  sometimes  raiseil  that  it  mijtht  he  inappropriate  to  apply  impedanet- 
concepts  ilevelopeil  for  TKM  transmission  lines  anil  hollow  waveguides  to  niicrostri|> 
lines  where  the  fields  no  longer  arc  constrained  to  a  finite  area.  This  question  has  heen 
dismissed  succinctly  hv  ( .etsin^cr'’  who  concludes  that  only  one  definition  of  microstrip 
characteristic  impedance  is  consistent  with  the  more  general  approach  based  on  wave- 
impedance.  which  •’ives  a  unique  result.  This  is  a  normal  definition  and  allows  appliea- 
tion  of  the  •'enerally  accepted  concepts  of  wave  impedance  and.  with  the  appropriate  ad¬ 
justments.  distributed  parameter  circuit  theory  . 

II.  TIIKORFTK.AI.  I)K\  KI.OPMF.NT 

If  the  line  has  finite  losses,  the  characteristic  inqieilance  and  propagation  constant.  7. 
can  he  calculated  from 


Zo  = 


n 


+  jwL 


+  jwC 


1 2a  I 


7  =  a  +  j/3  =  V(R  +  jwl.l  |(i  -F  jWC.l 


(21.1 


M.  A.  Wheeler,  “Transmission  Line  Properties  of  Parallel  Strips  Separated  by  a  Dielectric  Sheet."  II Tl  Trans, 
on  Microwave  Theory  and  Techniques,  Vol.  MTT-1  3,  No.  2,  pp  I  72-185,  Mar  65. 

T.  (*•  Bryant  and  J.  A.  Weiss.  "Parameters  of  Mierostrip  Transmission  l  ines  and  of  Coupled  Pairs  of  Micro- 
strip  Lines."  II  1  I  Trans,  on  Microwave  Theory  and  Techniques,  Vol.  MTT-1 6.  No.  12,  pp  1021-1027.  Dec  68. 

W.  .1.  Cictsinper.  “Mierostrip  Characteristic  Impedance,"  II  IT  Trans,  on  Microwave  Theor\  and  Techniques. 
Vol.  MIT-27.  No.  4,  Apr  79. 


when*  K  uml  (i  are  the  series  resistance  and  shunt  conductance  per  unit  length  of  the  lint 
and  a  and  0  are  the  attenuation  and  phase  constant  of  the  line.  With  the  appropriati 
modifications  to  account  for  normal  low-loss  transmission  line  and  other  reasonahh 
assumptions  given  hv  Brooke  et  al..7  an  alternate  equation  for  characteristic  inipcdanct 
results  in  the  form 

Zo=  el.  (.'5 

where  I.  is  the  inductance  per  unit  length  and  is  direetlv  affected  In  Both  the  resistance  ol 
the  line  and  frequency.  The  geometry  to  he  solved  is  shown  in  Figure  I. 


Figurel.  Transverse  geometry. 

For  generality  the  two  conducting  tapes  are  allowed  to  have  differing  widths:  however, 
they  will  he  maintained  parallel  and  with  their  midpoints  defining  a  plane  normal  to 
both.  This  is  a  simplification,  and  is  not  required,  hut  will  shorten  computation  time 
substantially.  This  configuration  can  he  used  to  represent  microstrip  examples  given  hv 
Vi  heeler’  and  Bryant  and  Reiss’  and  direct  comparisons  of  results  made  for  the  limiting 
eases  of  high  frequency  and  no  loss.  In  addition,  this  arrangement  permits  the  calculation 
of  the  parameters  of  antenna  feed  lines  where  the  widths  are  the  same  as  well  as  the 
caleulation  of  microwave  components  where  the  ground  plane  has  a  known  finite  width. 

Mathematically  subdividing  the  conductors  into  smaller  parallel  sections  is  accom¬ 
plished  as  show  n  in  Figure  2.  This  method  of  subdivision  is  arbitrary  and  is  retained  for 
consistency.7 


The  two  conducting  tapes  have  now  been  replaced  hv  In  thin  parallel  tapes,  each  of 
which  may  carry  a  different  current.  An  equivalent  circuit  of  the  transmission  line  then 
IfMiks  like  that  in  Figure  3. 

^  11  A.  Wheeler,  “Transmission  Line  Properties  of  Parallel  Strips  Separated  by  a  Dielectric  Sheet."  11  1!  Trans, 
on  Microwave  Theory  and  Techniques,  Vol.  MTT-1 3,  No.  2.  pp  1 72-185.  Mar  65. 
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T.  G.  Bryant  and  J.  A.  Weiss,  “Parameters  of  Microstrip  Transmission  Lines  and  of  Coupled  Pairs  of  Microstrip 
Lines,”  II* FI*  Trans,  on  Microwave  Theory  and  Techniques,  Vol.  MTT-1 6.  No.  12,  pp  1021-1027,  Dec  68. 

R.  L.  Brooke,  C.  A.  Uocr,  and  C.  II.  Love,  “Inductance  and  Characteristic  Impedance  of  a  Strip-t ransmission 
Line.”  journal  of  Research  of  NBS,  C.  Engineering  and  Instrumentation.  Vo!  7 1C.  No.  1.  Jan-Mar  67. 


3n+2 )  3n+l 


L  ■  -  — * 


_ _ ig  subdivisions. 


Mll  'l 
TOT— 1 


The  width  of  each  subsection  will  he  chosen  sufficiently  small  to  consider  the  current 
density  in  each  to  he  uniform.  Appendix  A  contains  the.  inductance  equations  which  can 
he  used  to  calculate  the  mutual  inductance  between  any  two  subsections  and  the  self  in¬ 
ductance  of  eaeli.  The  resistance  of  each  section  w  ill  he  the  tl.c.  resistance  calculated  from 
the  input  parameters  of  hulk  resistivity  anti  incremental  area.  The  incremental  area  is 
defined  b\  the  smaller  of  the  actual  tape  thickness  or  an  arbitrary  multiple  (an  input 
parameter)  of  the  skin  depth  anti  the  subsection  width.  In  this  paper,  the  upper  anti  lower 
tapes  are  each  divided  into  2n  etpial  width  subsections.  Brooke  et  al..7  8  discuss  alter¬ 
native  methods  used  to  speed  convergence  for  different  geometries.  Physical  symmetry 
produces  the  follow  ing  relations  between  the  currents  in  the  different  sections.  For  the  top 
tape. 

~  L>n  +  j.  j  =  I.2....11. 

For  the  lower  tape. 

F  —  I 

4  it  +  j  *  .'in  +  j.  j  =  1.2... .11. 

Tb<*  total  voltage  drop  per  unit  length  along  any  one  section  may  he  written. 

4n 

vk  =  rk  *k  +  j"  Mki>  •  {i) 

i=i 

Then*  are  4n  equations  which  may  he  written  in  matrix  for  as 


As  a  result  of  the  symmetry  conditions,  only  2n  of  the  equations  are  independent.  The  two 
groups  k  =  I  to  n  and  2n  +  1  to  3ii  are  equivalent  as  are  the  groups  k  =  n  +  1  to  2n 
and  Hn  +  1  to  4n.  The  order  of  the  matrix  equation  can  he  reduced  by  a  factor  of  two  hut 
it  should  he  noted  that  all  subsection  currents  affect  all  others  and  their  contribution 
must  he  reflected  in  the  final  equations.  Choosing  the  left  half  of  each  tape  as  the  in¬ 
dependent  equations  to  he  solved,  rewrite  (;>)  as 

7  R.  L.  Brooke.  C.  A.  Hoer,  and  C.  H.  Love,  "Inductance  and  Characteristic  Impedance  of  a  Strip-transmission 
Line,"  Journal  of  Research  of  NBS,  Inc.,  C.  Engineering  and  Instrumentation.  Vol  7 1C.  Jan-Mar  67. 

**  R.  E.  Brooke  and  J.  E.  Cruz,  “Current  Distribution  and  Impedance  of  Lossless  Conductor  Systems,"  IEEE 
Trans,  on  Microwave  Theory  and  Techniques,  Vol.  MTT-15,  No.  6.  Jun  67. 


Th<‘  following  definition*  arc  u*cd  In  >iinplifv  1 1 »<•  manipulation*: 

u>  =  K'  ||l)  +  io)MR-Vf 
<t>  =  o)M  K  1 
7  =  1 1,  -  lj„l 

6  =  ID . On.  I  . |  |l.>| 

\»  radial  mu  i*  not  allowed  and  the  line  i*  driven  a?  a  Lecher  “o-and-relurn  eireuil.  the 

total  circuit  in  the  one  (ape  tnu*t  cipial  the  total  current  in  the  other.  Thi*  cive* 

n  ’n 

L  \  *  £  \  e  0 

t  n  + 1 

n  2n 

E  >\ +  £  \ s  0 

I  n  +  1 

or  in  the  niutrie  notation.  U'inj;  1 1  d).  ( I  I).  and  (1.7) 

7  \  =  7to|K  +  4>V |  =  o  ilO) 

-yH  =  7to|<^>K  —  I)  so.  ,17, 

\-  all  the  currents  are  longitudinal  there  cannot  he  anv  lran*\er*c  voltage  difference'  on 
the  conducting  tape*  and  the  *ealar  voltage*  ohev  the  followin':: 

<i  =  =  1.2... .n 

•\ 

n  +  l....2n 


The  voltage  drop  around  the  loop  (down  one  tape  and  hack  the  other)  i*: 


Since  I  lie  value  of  the  voltage  drop  is  arbitrary  and  can  he  adjusted  to  any  quantity  hv 
ehaneine  the  applied  signal,  we  define  it  as  I  +  jO.  This  condition  applied  to  the  above 
results  produces 


The  column  matrices  K  and  K  can  now  he  written  in  terms  of  known  matrices  and  one 
scalar  uiikuow  n  each: 


K  =  —  61  +  e^1 

(Id) 

K  =  f|7' 

where  the  superscript  T  denotes  the  transposed  matrix,  (lomhiniii"  ( Id)  with  ( 16)  and  ( 17) 
and  sol  vine  for  e  and  f  produces 


<7|/'07T  Hy'l'QO '  )  +  (y\pOr  ) 

(71 py)2  +  (y\p<pyf  )2 

( 19) 

y\pOT  Ciyfty1 

y\p<f>yJ 

CO) 

Kvery  term  on  the  ri^lit  is  now  a  calculable  scalar.  The  results  obtained  can  he 
substituted  into  (Id).  ( Id),  and  ( 14).  and  the  current  in  each  section  determined.  Solution 
of)  Id)  and  (20)  is  performed  in  subroutine  (T  RR  in  a  straight  forward  fashion.  The  solu¬ 
tion  to  I  Id)  and  |  I  f)  mav  be  more  obvious  in  the  code  if  we  rewrite  those  equations  usin*r 
tin*  definitions  in  (15) 


V 

A  =  (K  +  </>K) 

(Ida) 

B  =  «*>K  -  K). 

(14a) 

Willi  the  currents  in  each  section  known,  the  approximate  inductance  and  resistance  per 
unit  length  can  be  calculated. 


Til.-  question  mijiht  lie  raised  as  to  whether  the  frequency  dependence  is  the  result  of  the 
finite  resistance  and  if  the  system  were  lossless  would  the  solution  for  la*ff  he  independent 
of  frequency.  An  examination  of  reference  (8)  shows  this  question  to  have  heen  addressed 
clearly.  For  clarity  and  completeness  the  argument  is  repeated  here. 


=  ]<•£>,  i,- 1.  =  -a-. 

which  leads  to  a  simple  form  for  the  current  matrix 

B  =  -L  VI"'  V. 

OJ 


l2.il 


(24| 


For  the  lossless  ease,  the  effective  inductance  is  found  to  he 


(25 ) 


which  still  appears  to  he  function  of  frequency. 


In  term-  til'  matrix  operation-  we  have 


T\  =  OB  =  1_0  \|-'  V. 
'  *  w 


t  —ft) 


ll ut-  resulting  in 


i  ..it  =  i/mr1  v  (27 1 

showing  t In*  Trt-i |i 1 1- 1 1<- x  independence  of  effective  inductance  umlrr  condition-. 

Numerical  results  obtained  through  a  programmed.  convergent  a|>|)n>\imali<»ii  for  I  .ell 
should  therefore  agree  with  capacitivilv  derived  re-lilt-  for  either  f—  oo  or  K  — *  o. 

111.  <:o\ykr<;k\<:k  and  <:i  rrknt  disthmu  tion 

The  com erjjenee  of  this  numerical  technique  i«  <le|iemlent  on  mam  factor-  including 
geometry.  frequency.  resistance,  ami  the  method  used  for  conductor  -ulidiv  i.-ioti.  Brookl¬ 
et  al.7  8  eni|doy  methods  of  tape  subdivision  using  uneipial  widths  to  speed  up  emner- 
lienee  and  provide  aeenraeies  better  than  0.1  percent  for  calculations  of  precision 
standard  Geometries.  To  provide  for  simple  adaptability  to  a  vurielv  of  Geometries,  the 
program  developed  to  calculate  mierostrip  properties  for  this  report  (  Appendix  B|  used 
only  equal  width  subdivisions.  Where  the  two  conductors  were  of  uneipial  width,  thev 
were  both  subdivided  into  an  equal  number  of  subelements  without  am  complex  scheme 
for  matching  the  current  distribution.  Figure  f  plots  the  normali/.ed  transverse  current 
densities  for  a  typical  mierostrip  with  equal  conductor  widths,  a  frequence  of  III7  II/.. 
conductor  resistivity  of  1.724x10"*  ohm-meters  and  an  impedance  of  about  .78  ohms.  Dif¬ 
ferent  plots  are  shown  as  the  number  of  subdivisions  are  increased  from  f  to  40  on  each 
conductor.  For  this  particular  ease,  convergence  (of  impedance)  of  better  than  0.5  percent 
was  achieved  for  u  =  1  I  or  22  subdivisions.  No  General  conclusions  should  be  drawn  from 
this  for  other  cases  or  Geometries.  Much  more  rapid  convergence  could  be  achieved  hv 
tailoring  the  suhdiv isioti  widths  to  each  particular  problem.  A  most  elegant  programming 
solution  would  be  to  adjus*  the  widths  through  an  iterative  process  of  equalizing  the  dif¬ 
ferences  between  adjacent  siihelcments.  ll  i-  not  obvious  that  such  a  scheme  would  result 
in  more  efficient  code  as  redundancy  in  impedance  calculations  would  decrease. 

If  the  conductors  arc  allowed  to  have  different  widths,  the  current  distribution  will 
show  a  markedly  different  pattern.  Figure  5  shows  three  overlaid  plots  for  one  ease  where 
the  wider  conductor  is  2.  .4.  and  5  times  the  width  of  the  other.  The  area  under  the  curve- 

^  R.  L.  Brooke,  (\  A.  Iloor.  and  (\  II,  Love,  “Inductance  and  Characteristic  Impedance  of  a  Strip-transmission 

Line,”  Journal  of  Research  of  NBS.  C.  Lnpineering  and  Instrumentation,  Vol  7 1C,  No.  I,  Jan-Mar  67. 

K.  I-.  Brooke  and  J.  I  .  Cruz.  “Current  Distribution  and  Impcdanec  of  Lossless  Conduetur  Systems.”  II  I  I 

I  rans.  on  Microwave  Theory  and  Techniques,  Vol.  MTT-15.  No.  6.  Jun  67. 
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tor  the  two  conductor*  should  he  the  same  for  each  case,  to  reflect  the  requirement  that 
the  total  current  he  the  same.  As  i*  expected,  the  current  on  the  wider  tape  crowds  to  the 
center  under  the  smaller  tape  and  reflects  the  partially  self-shielding  characteristic  of 
microstrip.  The  characteristic  impedance  obtained  for  these  three  cases  was  96.87  ohms 
(2/ 1 1.  92.21  ohms  (8/ 1 1.  and  89.69  ohm*  (.7/ 1 1.  (  hir  program  allow*  an  examination  of  the 
effect  ol  the  finite  w  idth  of  the  ground  plane  on  the  properties  of  microstrip  and  reveals  a 
surprising  dependence  even  for  ratios  of  1 0/1.  Impedance  value*  are  found  to  change  a* 
much  as  0.5  percent  for  a  10-pcreent  change  in  ••round  plane  width  at  a  nominal  10/1 
ratio.  \  recent  report''  goes  to  great  effort  to  establish  convergence  better  than  0.1  percent 
and  vet  is  based  on  an  assumption  of  infinite  width  for  the  ground  plane. 

IV.  COMPARISON  OK  RKSU.TS 

Two  cases  will  be  calculated  and  compared  with  results  produced  bv  other  author*.  An 
equal  width  ease  can  be  compared  directly  with  Wheeler's3  results,  and  an  unequal  width 
case  can  be  compared  with  the  results  of  Rrvant  and  Weiss'*  providing  the  larger  tape  i*  at 
least  10  times  the  width  of  the  smaller,  as  the  latter  reference  assumes  an  infinite  ground 
plane.  Only  the  high-frequency,  or  lossless,  case  will  be  considered  since  this  is  also  an 
assumption  of  the  references.  The  comparative  values  obtained  from  the  reference*  re¬ 
quired  interpolating  published  response  curves.  The  excellent  agreement  with  Wheeler3  i* 
more  than  would  have  been  expected  (Table  1).  The  slightly  high  bias  of  the  results  in 
Table  2  are  the  result  of  approximating  an  infinite  plane  with  a  10:1  width  ratio.  The 
program  calculations  have  been  rounded  to  the  nearest  ohm. 

The  current  distributions  produced  in  solving  for  the  results  of  Table  1  are  shown 
superimposed  in  Figure  6  for  values  of  w/h  from  I  to  80.  It  i*  clear  that  widening  the  two 
conductors  results  in  a  more  uniform  distribution  of  current  in  the  transverse  plane  and 
better  shielding.  This  i*.  of  course,  w  hat  one  should  expect  and  i*  used  by  Wheeler1  a*  the 
basis  for  a  wide  tape  approximation  and  to  establish  a  limit  for  the  effective  dielectric 

constant  of  K  and  propagation  velocity  of  v K. 


II.  A.  Wheeler.  ‘Transmission  l  ine  Properties  of  Parallel  Strips  Separated  by  j  Dielectric  Sheet.”  II  I  I  Trans, 
on  Microwave  Theory  and  Techniques,  Vol.  M  I  T-13,  No.  2,  pp  172-185,  Mar  65. 

T.  (i.  Bryant  and  i.  A.  Weiss.  “Parameters  of  Microstrip  Transmission  Lines  and  of  Coupled  Pairs  of  Microstrip 
Lines,”  II  1  1  Trans,  on  Microwave  Theory  and  Techniques.  Vol.  M  IT-16.  No.  12.  pp  1021-1027.  Dec  68. 

R.  11.  Jansen.  “High-speed  Computation  of  Single  and  Coupled  Microstrip  Parameters  Including  Dispersions. 
High-order  Nodes.  Loss,  and  finite  Strip  Width,”  II  I  f  Trans,  on  Microwave  Theory  and  Techniques.  Vol. 
MTT-26,  No.  2,  I  eh  78. 
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Table  1 .  Comparison  of  Results  for  hqual  Width  (K=l ) 


w/h 

Wheeler3 

This  Method 

0.3 

315 

313 

0.4 

279 

280 

0.5 

252 

254 

0.6 

232 

234 

0.7 

216 

217 

0.8 

202 

202 

0.9 

189 

190 

1.0 

178 

179 

3.0 

87 

87 

10.0 

32 

32.5 

^  11.  A.  Wheeler,  “Transmission  Line  Properties  of  Parallel  Strips  Separated  by  a  Dielectric  Sheet,”  ILLL  Trans,  on 
Microwave  Theory  and  Techniques.  Vol.  MTT-  1 3,  No.  2.  pp  172-185,  Mar  65. 


Tabic  2.  Comparison  of  Results  for  Unequal  Width  (K=l ) 


w,  /h 

W 2 / W ,  =  10 

Bryant  &  Weiss4 

This  Method 

0.6 

156 

160 

0.8 

140 

142 

1.0 

127 

128 

2.0 

87 

90 

3.0 

(.7 

70 

4  I .  (1.  Hr >  Jilt  ;mc!  .1.  A.  Weiss.  “Parameters  of  Microstrip  I'ransmission  Lines  and  of  Coupled  pairs  of  Microstrip 
Lines.”  II  1  1  Trans,  on  Microwave  l  heorv  and  Techniques.  Vol.  M  l  1-16,  No.  12,  pp  1021-1027.  Dec  68. 
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V.  IWKAWKTKK  FHK(,)1  KNCY  DKPKNDKNCK 


Tilt-  transverse  current  distribution  is  a  flint-lion  of  frequency  anti  lilt-  other  physical 
variables.  The  effective  intliielanee  per  unit  length  anil  allenuation  tlerivetl  from  the  eur- 
renl  ilislrihiition  will  also  he  fretpienev  tlepentlenl  anil  fully  preiiietahle  with  this  pro¬ 
gram.  Tilt*  alleniiatioii  tlerivetl  ht-re  is  only  the  ohmic  loss  in  the  eonihiclors  anti  tines  not 
account  for  loss  anti  tlispersion  in  (he  tlieleelrie  suhslrale. 

To  examine  the  behavior  of  both  etpial  anti  unetpial  tape  microstrips  as  a  function  of 
fretpienev  we  select  a  w/h  ratio  of  5/1  for  the  etpial  case  ami  w(/h  ratio  of  4/1  for  the 
unetpial  case,  with  w./w^  =  10  to  assure  re asonahlt-  similarity  to  typical  components 
usetl  hv  iutluslrv  in  tin-  manufacture  of  microwave  components  anti  antennas. 

The  frequency  tlepentleiice  of  the  transverse  current  distribution  for  these  two  cases  pro¬ 
vides  insight  into  w  here  the  transition  occurs  anil  the  cause  for  ohmic  losses  w  hich  can  In¬ 
substantial.  Figure  7  is  a  parametric  presentation  of  the  transverse  current  distribution 
for  an  equal  conductor  width  mierostrip.  (w  ,  =  w\,  =  .01  m.  h  =  .002m:  p  =  1.7241  x 
10'*)  The  major  frequency  dependence  occurs  between  1()4  Hz  anti  10s  Hz.  At  a  frequency 
of  5  x  10s  Hz,  the  inductance  anti  characteristic  impedance  have  stabilized  although  the 
ohmic  attenuation  w  ill  continue  to  increase.  Figure  8  shows  the  frequency  response  of  the 
unetpial  ease  which  shows  a  much  greater  tlepentleiice  but  is  still  sett  let  I  down  by  10*  Hz 
(w,  =  .01m,  w2  =  0.1m,  h  =  .0025)  (note  that  the  display  program  plots  distribution  at 
full  frame  size  when  w\,/w(  >  10). 

If  we  examine  the  frequency  dejiendence  of  these  two  types  of  geometries  we  can  gain 
an  understanding  of  their  behavior  anti  the  effect  of  geometry  on  it.  In  Figure  9.  we  plot 
the  characteristic  impedance  of  an  etpial  width  mierostrip  w  ith  a  dielectric  of  1  as  a  func¬ 
tion  of  frequency.  There  is  a  modest  decrease  from  the  d.e.  value  (60.4  ohms)  to  the  high 
frequency  limit  (58.55  ohms).  The  change  occurs  in  the  fretpienev  range  of  5  x  UP  to  5  x 
10s  Hz.  If  the  matierial  used  were  a  better  conductor,  this  shift  region  would  be  translated 
to  the  left.  Ill  sharp  contrast  to  this  is  the  behavior  of  unetpial  width  mierostrip  where  a 
much  larger  change  is  observed  with  increasing  fretpienev.  Figure  10  is  a  typical  plot  of 
one  ease  where  the  low  frequency  impedance  of  124.3  ohms  is  transformed  to  a  high- 
frequency  impedance  of  57.33  ohms.  The  transition  commences  at  a  lower  frequency  and 
is  essentially  complete  at  a  frequency  of  10s  Hz.  The  more  ratlit-al  dependence  on  fre¬ 
quency  is  the  result  of  large  shifts  ill  ground  plane  currents  to  concentrate  under  the 
smaller  conductor. 

Figure  1 1  plots  the  characteristic  impedance  of  etpial  width  mierostrip  for  four  dif¬ 
ferent  effective  wavelengths.  X  .  X  12.  X  /3.  anti  X  /4.  as  a  function  of  the  ratio  of  conthic- 
tor  width  to  separation.  Wheeler’s  wttrk  can  lie  usetl  to  estimate  the  relative  fill  factor,  q. 
which  can  be  usetl  to  determine  the  effective  dielectric  constant  K'.  Following  W  heeler  we 
have. 
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Figure  7.  Frequency  dependence  of  current  distribution  for  equal  width  microstrip. 
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TRANSVERSE  CUT 

Figure  8.  Frequency  dependence  of  current  distribution  for  unequal  width  microstrip. 
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U1I0’  =  cM=lM 


Figure  9.  Equal  conductor  microstriop  characteristic  impedance  as  a  function  of  frequency. 


.0025 


FREQUENCY  (Hz) 

Unequal  conductor  microstrip  characteristic  impedance  as  a  function  of  frequency. 


K  =  I  +  q  <k-l|.  (2‘>) 

w here  K  is  lilt-  nominal  dielectric  constant  of  (In-  dielectric  sheet  between  the  conductors. 
I  sin«  Wheeler's  Fi«un*  l>  yields  an  approximate  value  for  <|.  The  effeetive  wavelcneth 
can  lie  reasonably  a|>|iro\imaIed  In  equation  (30). 


(3<»| 


An  assumption  used  here  is  that  the  limit  of  K'  lor  verv  narrow  tapes  is  |k+  1 1/2  or  an 
average  w  ith  free  spaee  and  the  limit  for  verv  w  ide  tapes  is  just  K.  the  dieleetrie  eonstant 
of  the  material  between  the  eonduelors.  Figure  12  presents  the  results  for  the  ease  of 
unequal  width  mierostrip. 


CHARACTERISTIC  IMPEDANCE  Of  UNEQUAL  WIDTH  MICROSTRIP 
AS  A  FUNCTION  OF  AND  EFFECTIVE  WAVELENGTH 


Figure  T2.  Characteristic  impedance  of  unequal  width  microstrip  as  a  function  of  W|  /h  and 
effective  wavelength. 


VII.  CONCl.lMON 


111  this  report  wr  have  shown  the  theoretical  development  of  a  method  lo  calculate  the 
inductance  | M-r  unit  length  and  ohmic  losses  of  micro-trip  lint-  with  general  t  rio-  sectional 
geometry.  A  program  dew-loped  lo  apply  thi-  teeliniipie  wa-  exercised  for  simple  ca»es 
and  lilt-  results  found  to  agree  with  those  of  other  authors  using  eapaeitanee-hased  solu¬ 
tions.  This  approaeh  is  unique  in  that  it  dircctlv  provides  tin-  transverse  eurrenl  distrihu- 
( ion  and  ohniie  attenuation  at  all  frequencies.  It  can  provide  new  insights  into  faelors 
which  cause  loss  ill  microwave  eonipoiients  and  help  explain  the  effective  behavior  of 
currents  oil  extended  antenna  structures.  A  clear  understanding  of  antenna  feedlines  and 
radiatin'!  elements  can  he  reached  onlv  if  the  current  distribution  is  known. 

This  particular  approaeh  has  a  special  advantage  in  the  ease  with  which  it  can  he 
extended  to  other  geometries.  It  must  he  cautioned  that,  although  the  theory  is  straightfor- 
ward.  there  are  a  tuple  opportunities  to  lose  the  solution  in  writing!  the  program. 

This  report  was  restricted  to  the  odd  mode  id'  single  element  inierostrip  heeause  the 
specific  application  for  which  it  was  developed  required  onlv  this.  It  could  he  extended  to 
even  mode  and  coupling  calculations  if  the  total  currents  can  he  defined  in  sufficient 
exactness  to  allow  solution  of  the  matrix  equations.  The  results  already  shown  here  can 
have  utility  in  producin'!  inierostrip  w  ith  minimal  ohniie  losses  and  reduced  cost  bv  plac¬ 
in'!  low  resistance  materials  (cold  or  silver)  only  in  the  area  id'  the  conductors  with  high- 
current  densities.  In  the  ease  of  unequal  width  inierostrip  this  means  a  few  millimeters  on 
the  outer  edges  of  only  the  smaller  conductor. 

VII.  PROGRAM  \OTKS 

This  method  has  been  presented  previously.7  hut  to  our  knowledge  has  not  been  used  in 
the  Id  years  since  publication.  Perhaps  the  more  universal  availability  of  large  high¬ 
speed  computers  and  the  increasing  acceptance  of  numerical  solutions  will  provide  en¬ 
couragement  for  its  use.  It  is  in  this  spirit  that  the  program  is  provided.  Appendix  (’.  con¬ 
tains  sample  runs  to  assist  in  program  checkout. 

The  program  runs  interactively  through  Intercom  on  a  ('.DC.  6600  using  Scope  3.1.  The 
code  provided  makes  a  modest  attempt  to  minimize  core  required  and  lo  reduce  run  time, 
hut  we  eschew  heroic  efforts  which  make  modification  and  transportability  difficult. 

Subprograms  Not  Provided: 

The  follow  ing  subprograms  are  called  from  subroutine  C.l’RR  and  are  found  in  Library 
3  of  the  IMSL  Fifth  edition.  November  1075: 

7  H  !..  Brooke.  C.  A.  Ilocr.  and  ('.  It.  love,  “Inductance  and  Characteristic  Impedance  of  a  Strip-transmission  lane." 

Journal  of  Research  of  NBS.C.  I  npneerim- and  Instrumentation.  Vol  7H  .  No.  1.  Jan-Mar67. 
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\  Ml  LFH  Matrix  Multiplication  (full  bv  band). 

\  Ml  I.FF  Matrix  Multiplication  (full  bv  full). 

I.IW2F  Matrix  Inversion. 

I.F(,)T2F  Solution  of  Linear  Filiations. 

Subroutine  I’l.OTXY  is  a  inaehine  anil  terminal  specific  plotting  subroutine  ealled 
from  subroutine  1*1.011.  It  is  not  necessary  to  inelinle  IM.OTW.  or  IM.OTI.  to  use  the 
program.  The  eonunents  whieh  proxiile  instructions  for  using  I’LOTXN  are  inebnleil  so 
that  the  user  max  interpret  the  calls  in  a  manner  compatible  with  his  nun  graphic 
terminal  if  desired. 

Function  Second  —  See  notes  below: 

Subroutine  ( ion  nee  (file  name)  associates  the  named  file  with  the  terminal  device  and 
max  be  an  unneeessarx  function  in  some  installations. 

System  Specific  Subprogram: 

The  following  subprograms  are  specific  to  the  current  ('.DC.  operating  system  in  use  at 
MFRVDC.OM  and  will  normailx  require  adaptation  or  removal. 

Function  Second  —  This  function  returns  the  central  processor  lime  consumed  by  the 
job  since  start. 

Subroutine  C.PIITIME  —  This  is  a  compass  subprogram  called  from  Function 
TIM  FI..  It  returns  the  time  remaining  before  the  system  will  abort  the  job. 

These  subprograms  are  used  to  provide  the  interactive  user  with  information  which  max 
be  used  to  manage  the  consumption  of  resources.  Vi  here  resources  are  not  a  constraint, 
dummy  subroutines  may  be  provided  or  the  calling  statements  removed.  Retention  of 
these  functions  is  generally  recommended  for  both  use  and  program  modification  as  the 
potential  for  resource  consumption  is  significant.  The  use  of  List-Directed  read  and  write 
(print)  statements  such  as  RF\D  *.  IN  arc  freely  intermixed  with  formatted  l(>.  If  no 
equivalent  feature  is  provided  bv  the  user’s  sVstem.  formatted  statements  will  have  to  be 
substituted.  Many  such  -tatements  which  were  inserted  during  program  development  have 
been  retained  as  comments. 


File  l^sage.  The  principle  output  lilt-  is  <)l  TPl  'I',  which  is  equivalent  to  TAPK6.  ami 
is  the  terminal  for  interactive  users. 


The  only  input  file  is  INPl  T.  which  is  the  terminal  for  interactive  users.  TAPK2  is  an 
alternate  output  file  used  to  save  specified  results  for  later  printing  under  interactive  con¬ 
trol  (not  the  terminal). 

Program  l  se.  The  program  is  uearlv  self-e\planatorv  when  run  iuteraelivelv.  The 
user  enters  at  least  the  first  two  characters  of  a  key  word.  If  the  kev  word  represents  a 
variable,  the  value  will  he  requested.  If  the  kev  word  represents  a  command.  (I)isplav 
variables.  Print  results.  (»0  calculate.  S  TOP,  instructions.  SAY  K).  that  function  is  per¬ 
formed.  Samples  of  operation  are  included  (  Appendix  (!).  Vi  hen  parametric  design  studies 
are  made,  the  user  would  be  wise  to  conserve  resources  hv  being  sensitive  to  which 
variables  require  which  calculations.  Changes  of  basic  geometry  require  eompletelv  new 
solutions,  while  frequency,  tape  thickness  and  resistivity  do  not  require  recalculation  of 
the  mutual  inductance  matrix.  Changes  of  dielectric  constant,  power  factor,  and  com¬ 
mand  inputs  are  trivial.  The  user  is  prompted  hv  messages  indicating  (d*l  time 
consumed. 

Specific  Subprogram  Notes: 

Program  STHIPZM  includes  only  trivial  ealeiilations  and  normal  output,  but  not  the 
input  functions.  A  somewhat  convoluted  method  of  specifying  arrays  is  introduced.  This 
method  is  helpful  during  program  development  and  modification  as  it  permits  flexible 
creation  of  arrays  and  core  management  without  substantial  changes  of  code,  but  at  ex 
edition  time  expense.  Two  dimension  arrays  of  dimension  (\1*\.  \1*\)  are  stored  in 
order  in  blank  common.  \  is  the  input  variable  spccifving  the  number  of  subdivisions  in 
a  half-tape.  \l  =  I.  2.  f.  H.  12.  If>.  The  number  of  arrays  of  (M*\.  \l*\)  dimension  is 
stored  in  vector  VA  where  the  index  to  VY  corresponds  to  the  M  implied  KXYYIPI.F: 
VY  =  2.  1.  0.  0.  .A.  0  means: 

2  arrays  of  Dimension  (\.\)  followed  hv 
I  array  of  Dimension  (2*V  2*\|  followed  hv 
■A  arrays  of  Dimension  |I2*V  I2*\) 

\  similar  scheme  for  the  storage  of  vectors  of  dimension  (\1*\ )  is  used  with  the  number  of 
vectors  of  dimension  (\1*N|  stored  ill  vector  NY  .  Y'cctor  storage  space  follows  the  arrav 
storage.  Finally,  the  particular  array  within  a  group  of  arravs  (or  vectors)  of  the  same 
dimension  may  lie  indirectly  specified.  Vector  ISA  is  used  to  perform  this  function. 
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Subroutine  IO  contains  the  input  code  and  may  be  replaced  by  those  who  prefer  only 
batch  processing.  The  logic  provided  attempts  to  reduce  execution  time  bv  requiring  only 
those  calculations  which  arc  necessary.  Subroutine  (T'RR  contains  the  principle  code 
which  achieves  the  purpose  of  this  paper.  Note  that  double  precision  is  used  in  the  sum¬ 
mations  required  for  the  solution  for  K  and  F.  This  subroutine  uses  core  beyond  the  work 
area  matrix  as  noted  in  the  program.  If  N  vector  space  is  required  to  survive  execution. 
storage  allocation  must  be  modified. 

Function  SZ1N1)  performs  the  mutual  inductance  calculation.  Execution  time  is  poor 
and  evaluation  requires  high  precision  as  the  difference  between  similar  values  occur. 
When  installed  on  a  machine  with  less  precision,  care  must  be  exercised  in  determining 
the  range  of  input  variables  to  be  permitted.  (It  is  possible  to  increase  N.  subdivisions, 
and  decrease  the  precision  of  the  result.)  Subroutine  V1I)X  takes  advantage  of  the  half 
plane  symmetry  and  the  equal  width  subdivision  of  tapes  to  permit  a  reduction  of  the 
number  of  mutual  inductance  calculations. 

Subroutine  MUT  calls  function  SZIM)  to  fill  a  work  space  with  the  unique  mutual 
calculations  required,  and  then  fills  the  complete  mutual  matrix  with  appropriate  values. 

Function  IN  A  provides  the  index  to  the  location  of  any  element  in  any  array  specified 
in  the  manner  described.  The  function  also  provides  an  entry  which  determines  the  max¬ 
imum  number  of  subdivisions  based  upon  the  amount  of  blank  common  available  and 
the  arravs  defined  by  vectors  N  V  and  NY.  These  values  are  passed  in  labled  common 
STOR. 


APPENDIX  A 


BASIC  EQUATIONS  FOR  INDUCTANCE  IMPEDANCE 


The  mutual  inductance  per  unit  length  between  two  long.  thin,  parallel  tapes 
such  as  shown  in  the  figure  below  can  be  obtained  from  equation  8  of  reference 
(10)  and  is 


A/. 


0.2 

ad 


ln(/)2  +.v2 )  -  xP  Tan  '_L 
P 


I:  -  a,  I ■’  +  d 
(x) 

E  +  </  -  a,  E 


+0.2|  In  21+  '/il.jull/m 


(Al) 


where  the  limits,  which  have  been  retained  for  compactness,  are  substituted  as  follows: 

s.s,  4 

\f(x)\ i-i  I*1#*,). 

s:s4  1=  I 

The  origin  coincides  with  the  left  edge  of  one  tape.  If  the  left  edge  of  the  second 
tape  is  in  any  quadrant  other  than  the  first  as  shown,  one  or  both  of  the  values  of 
E  and  P  will  be  negative.  The  self-inductance  per  unit  length  of  a  long  thin  tape  is 

Lt  =  0.2  In  J_  +0.2  |  In  2/  +  '  :| ,  /all/m.  (A2) 

a 


In  both  of  these  expressions  it  is  assumed  that  the  current  density  is  uniform  through 
out  the  conductors. 


to 


C.  A.  Ilocr  and  C.  II.  Love.  I.xact  inductance  equations  for  rectangular  conductors  with  applications  to  more  com* 
plicated  geometries,  J.  Res.  NBS69C2  (Lngr.  and  Instr.),  No.  2.  127-137  (Apr.- June  1965). 


Cross  sectional  view  of  two  long  parallel  thin  tapes. 


Note  that  the  term 


0.21  In  21  +  Zi]  (A3) 

appears  in  both  \1t  and  L{  as  the  only  term  involving  length.  However,  if  M  and  / 
are  substituted  into  (8).  the  terms  involving  length  will  exaetly  cancel.  Substituting 
(A  1 )  and  (A2)  into  (8)  gives 

in  8^  f 

vk  =  rk'k  +  +  '°-2(  1  n  +  ^))2J l,  ' A4) 

/=  i  /=  i 

where  the  M'kl  are  the  expressions  in  (Al)  and  (A2)  with  the  (A3)  term  excluded. 

But  the  last  term  is  zero  because 


/=  l 


That  is.  the  total  current  in  the  outer  conductors  is  equal  and  opposite  to  the  total 
current  in  the  inner  conductors.  Therefore,  the  Mk  in  (8)  may  be  calculated  from 
(Al )  and  (A2)  with  the  (A3)  term  excluded. 


APPENDIX  B 

PROGRAM  FOR  CALCULATION  OF  APPROXIMATE  CURRENT 
DISTRIBUTION,  INDUCTANCE,  AND  IMPEDANCE  OF  MICROSTRIP 


o  o  o 


PROGRAM  STRIPZMi TAPE2=65,  INPUT=  o5  ,OU TPUT=65/ 1 *♦  0  ,  T A PE6  = 

COUTPUT1 

C  THIS  PROGRAM  USES  THE  METHOD  PRESENTED  IV  INDUCTANCE  »ND  CHARACTERISE 
C  IMPEDANCE  DP  A  ST- IP  TRANSMISSION  LINE,  R.L.  3R0DKE  £1  Al,  JOURNAL  OF 
C  REStARCH  OF  THE  NATIONAL  BUREAU  OF  STANDARDS-  C. ENGINEERING  AND 
C  I NSfRUMENTATl ON  VOL «  7 1 C ,  NOl,  JAN-MAR  1967 

C 

C  MATRIX  INVERSION  AND  LINEAR  EQUATION  SOLUTION  SUBROUTINES  ARE  -OUND  IN 
C  LIBRARY  3  OF  THE  INTERNATIONAL  MATHEMATICAL  AVJ  STATISTICAL  LIBRARIES 
C  FIFTH  EDITION,  NOVEMBER  1975.  IOGT=0  MILL  PERFORM  ITERATIVE  IM^ROVEMEN 
C  T  OF  THE  SOLUTION  ANO  REPORT  THE  NUMBER  OF  DISETS  UNCHANGED  AFTER 
C  TM PRTv EHTNT.  IF  IOGT  not  ZERO  ON  INPUT,  IMPROVEMENT,  IF  REQUIRED,  IS 
C  ATTE-PTED  TO  I DGT  DIGETS. 

C  THE  PLOT  SUBROUTINES  ARE  MACHINE  SPECIFIC  AT  MERADCOM  AND  NOT  INCLUDED 

C 

C  THE  FUNCTION  INA  PERMITS  SPECIFICATION  OF  ANY  (I,J>  E-EHENT  IN  ANY  <KI 
C  SPECIFIED  ( N*M I  *  *  2  ARRAY  OR  N*M  VECTOR,  WHERE  M  HAY  BE  1,2, 3,12,16. 
C 

C  STORAGE  IS  IN  BLANK  COMMON,  ARRAY  A<ICDRE»  WITH  THE  SQUARE  ARRAYS 

c  first,  then  the  vectors. 

c 

C  THE  NUMBER  OF  EACH  ARRAY  TYPE  DEFINEO  IS  STORED  IN  NAIIH)  AND  VV(IM) 

C  WHERE  IM= 1 , 6  ( SPEC  I F I  NO  THE  M  IN  ORDER  GIVEN  ABDVE) 


C 

C 

complex  rpjhl,gpjwc,zcomp,gam 

COMMON  //A (27326) 

COMMON  /STOR/N,MAXN,ICORE,NV«6> ,NA(6>  ,MM(6> 
EXTEPNAL  INA.INV 


THE  ASSIGNMENT  OF  SPACE  FOR  ARRAYS  IS  MADE  THROUGH  THE  VECTOR  ISA  AND 
IS  AS  FOLLOWS 

ISA  IS  A  VECTOR  SPECIFING  THE  MATRIX  NUNBtR  FDR  VARI43LES  AS  FOLLOWS: 

M  MUTUAL  INDUCTANCE  MATRIX  IS  INAI2,ISA<1> 

PHI  MATRIX  INA (2  ,  ISA ( 2 ) 

PS  I (INVERTED)  MATRIX  INA(Z,ISA(3) 

PSI  MATRIX  INAC?,ISAU> 

WORK  AREA  MATRIX  INA(2,ISAC5* 


Q* *********************** ****************************** ****** 

C  THE  WORK  AREA  EXTENDS  6*N  BEYONO  I NA ( 2, I S A t 5 » , N2, V 2 1  AND  IS  PRDVIDEO 
C  BY  3  2  N  VECTORS  DEFINEO  BELOW,  IF  N  VECTORS  ARE  OErINED 

ELSEWHERE,  THEY  WILL  ef  WRITTEN  OVER,  THE  AREA  MAY  9E  USED 
FOR  VOLATILE  STORAGE, 


C  TWO  2*N  WORK  AREA  VECTORS 

C  THE  THIRD  ?  N  AREA  IS  USEO  FOR  1/R  INV(?,ISA(6) 

C 

C  A  VECTOR  INV<2,ISA(7» 

C  B  VECTOR  TNV(2,ISA(8) 

C  C  VECTOR  INV(2,ISA(9) 

C 

DIMENSION  ISAJ9) 

COMMON  /STRIPS/FPEQ*SEP»W12> , E W ( 2 ) , T HK.SKO^S , E THK 
C * RHO » EUC »PO W , I OG T  ,IER,PI2 
DATA  N/2/ 

DATA  ICQRE/2 78267 
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DATA  NA/0, 5,4- 3/ 

DATA  NV/O,  6,4*0/ 

OATA  MM/1,.2 ,4,  8,12 , 16/ 

DATA  IU,IUPR/2*6/ 

DATA  COND/C./ 

C 

DATA  ISA/1,2,3,.,,5,3,4,5,6/ 

CALL  CONNEC  (5L INPUT ) 

CALL  C0NNECI6L0UTPUT) 

PRINT  , 'ENTER  I  FOR  INSTRUCTIONS,  01  FOR  DEFAULT  VALUES  * 

TP  =  3 

CALL  ISETN <  N ) 

SOTO  ? 

2  ir=: 

3  CALL  IO(IR) 

C  RE TU°NS  AFTER  VARIABLES  ARE  READ  AND  A  COMMAND  INPUT  PROVIDED  IN 
C  SUBROUTINE  10.  VARIABLE  IK  IS  RETURNED  TO  CONTROL  THE  EXECUTION 
C  REQUIRED. 

GO. 0(13  0  ,1.1, 102, 1..3, 104, 1*5, 106), IRU 
C  NO  CHANGE  .REPEAT  BASIC  OUTPUT 
103  GOTO  1301 
C  DO  ONLY  LOCAL  CALCULATIONS 
101  GOTO  1032 
C  DO  CURRENTS  BUT  SKIP  MUTUAL 
132  I SK=  1 

GOTO  1333 
C  00  CURRENTS 
103  rs<=4 

1003  TIN=SECOND(CP) 

CALL  CURR!N, ISA, SUMA, SUMS, ISO 
TREQ=SECONO!CP)-TIN 
PRINT  910, TREQ 
NPM=TIMEL(TL)/TREQ 
IF(NRM.LE.8)PRINT  9 11 , TL , NRM 
9il  FORMAT!*  YOU  HAVE  *,F8.3,*  CPU  SECONDS  LEFT*/ 

C  *  ABOUT  ENOUGH  TIME  FOR  *,I3,*  MORE  RUNS  LIKE  THIS  */) 

910  FORMAT!*  THIS  RUN  REQUIRED  *,F6.2,*  SECONDS*/) 

GOTO  1.02 
C 

C  PRINT  CURRENTS 
105  LDCC=INV<2,  ISA<9)  ,1,11-1 
L0C=INV!2,ISA<7) ,1,1)-1 
LOCO = I NV! 2, ISA (8) ,1,1)-1 

WRITE (IU,9uQI  IA(LOC*I)  ,A!  L0C2+I)  ,«UOCC*I),A(lOCM*II, 

CA <L0C2*I»N) ,A  !L0CC  +  I»N1 ,I=1,N) 

903  FORMAT (T33 , *CURRENTS*/T1G*L0WER  TAPE*,T4C, 

C’UPPER  TAPt*/T3,*REAL-  ,T15,*-U  ( I  MAG ) ♦ , T 2 7 ,  "4  A G NI T UD E*  , 

C,T3  9,*REAL*,T51,*-U  (I  MAG) * , T63 , *  MAG N I TU DE * / 

C4C6(1X, 1PE11.4)/) ) 

TF! IU-IUPR) 1062,2 

C 

C  DO  SOME  TRIVIAL  CALCULATIONS 
C 

1002  CONTINUE 

HF=dI2*FREQ 

SAP=3UMA*SUMA+SUMB*SUMB 

CLEFF=SUMB/<MF*SAB) 


» # 
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of.  >*r.o  4.  r.oo 


ftEFF^SUMA/!  SA8> 

V  ELC =  2. 997925 ES/SQRT (EOC) 

CCEFFs  1./ (CLEFF* V£LC**2> 

CONO=HF*CC£F£' POM 
ZIMP- V£LC*CL£FF 

1031  W ft  I TE  < I U ,9b  3  > C  L  £ F  F , REFF ,Z  I MP.CCEFF,  ^UnD 

963  FORMAT ( *  EFFECTIVE  INOUCTANC£=  *,1PG11.4,*  HEM  E  R VS/M*/ 

C*  EFFECTIVE  RESISTANCES  »,1PG11.4,»  OHMS/M*/ 

C*  CHARACTERISTIC  IMPEDANCES  *,1PG11.4,*  OHMS*/ 

C*  EFFECTIVE  CAPACITANCES  *,1PG11.4,*  FARADS/M*/ 

C*  SHUNT  CONDUCTANCES  *,1PG11.4,*  MHDS/M*/ 

C)  '  ~  '  “ 

RPJHL=CMPLX (REFF,HF' CLEFF) 

GPJHC=CMPLX (COND,HF*CCEFF) 

PRINT*, 'MF, REFF, CLEFF, RPJNL.GPJWC-  *,HF, REFF, C LEF^ , RPJWL *G» JMC 
PRINT  +, *MF,VELC,Cl£FF,HF/(CLEFF*VE.C**2) • ,Hr , VELC , CLEFF  ,HF/ ( CLEFF 
C* VELC**2  ) 

7C0HP=C?(IRTT  RP  JHL/GPJMCT  —  — —  -  - 

GAM=CSQRT( RPJML’GPJMC) 

HRIT£(IU,972IZC0MP,GAM,8.686*REAL(GAM>  ,  3  60  .  *  A  I  MAG  C  GAM) /PI2 
72  FORMAT!*  COMPLEX  IMPEDANCES  * , IPG  13.  4 ,10 H  OHMS  J.1PG13.4, 

C5  H  OHMS/*  PROPAGATION  CONSTANTS  ♦  ,1P3 13. 4, 1 3 H  NEPERS/M  ,1PG13.4, 
C,10H  RAO  IANS/M  / 

CT77,  IPG  1374 , 1 3H  OB/H  riPGTlT.  47T3 FT1TE5 RETSTM  7 

C) 

IF(IU-IUPR) 1361,2 
PLOT  CURRENTS 

04  CALL  PL0TI(A(INV<2,ISA<9» ,  1 ,1 ) > , A (IMA (2 , 2, 1 , 1 1 »,2*N,M» 

~  — GOTO  "2  -  -  - - 

SAVE  ON  TAPE2 
106  IU=2 

GOTO  1301 

1061 _ GOTO  135  _ _ _ 

IU=IUPR 
IR  =  0 
GOTO  2 
ENO 


SUBROUTINE  10 ( I ft  > 


C 

C  RETURN  COOES  AS  FOLLOWS 
C  IR  =  0  NO  CHANGE 

C  1  ONLY  MAIN  PROG  CALCS  NECESSARY 

C  "  2  CURRENT  NECESSARY,  BUT  NOT  MUTUALS 

C  3  COMPLETE  CALCULATION  NECESSARY 

C  *♦  PLOT  CURRENTS 

C  5  PRINT  CURRENTS 

C 

EQUIVALENCE  (  F  REG  ,  V  < 1 )  ) 

COMMON  /STRIPS/FREQ, SEP, Wt2j ,EW(2) »THK,SKDPS, ETHK* RMO 
C,EOC.POW,IOGT, IER,°I2 
COMMON  /ST0R/N,IDUM(2G| 

Ly 

DIMENSION  W IN ( 2 ) 

DIMENSION  Q<21),KW(2i),UN«21»,V(12l 
C 

OA  T  A  NE/12/ 

DATA  Q/ 1  OH  ELEMENT  S  , 1 xHFPEQUFNGY  , 1 0 HSE3 AP ATI  ON  »6HHI DTH 1 , 

CbHWIDTH2,lJHELE.WlOTHl,lGHELE.WlDTH2,loHTHIC<NESS  ,10MSKIN  DPTHS 
CtlCHEFP. THICKS, 3HPH0,1GH0IE.CCNS  , 

C12HP0WER  FACT  , 10 H PESO LUT ION  , 2HG0, 4  H  PLO  T , 

CD  HFRINT ,4HST0P,1CH  DISPLAY  , I UH INSTRUCT  YITHSAVE  TAPE?  7 
DATA  KW/1HN,2HFR,2HSE,2HWI,2HWI,1H  ,1H  , 2HTH , 2H SK, 1 H  , 2hRH , 2 HED , 2H 
Cp0,2HI0, 2HG0,2HPL , 2hPR, 2HS T , 2H0 I , 1HI , 2HSA/ 

DATA  UN/1h  ,2HHZ,6- (6HMETERS)  ,  1H  ,6H  M  ETE  RS ,  1 G  H  OHMS  *9  ET  ER 
C.IOHRELATIVE  ,5h(SIN),6HDIGETS,7*(1H  )/ 

DATA  V  /  1  u  a  •  ,.01,.«1,.01,u,,0.,.U'jj25,1.,o., 

-  TI.TT9IE-  38,  1 .  ,  0,  / 

DATA  IDG T, IER,PI2/- ,0 ,6.283185339/ 

DATA  Ip/0/ 
r 

vy 

5  PRINT  *, 'ENTER  KEYWORO  ' 

PE  A  0  9U,IA 

'901  FORMAT  f  A  21  - -  - 

L=Nt*9 
DO  1 j  1=1, L 

I F (  IA.EQ.KWdl  (GOTO  15 
10  CONTINUE 

GOTO  5 

15  TFtT.CQ Y6TG0TO  5  -  .  .  -- - - 

IF(  I.  EQ.DGOTO  28 

IF( I.GT. l*NE(GOTO  25 
IF ( I. EQ.4) GOTO  20 
PRINT  ' , 'ENTER  VALUE  ' 

READ  *,V(I-1» 

IFd.EQ.  3TIR  =  3  - - 

IF< I.LT. 12) GOTO  17 
IF(Ift.LE.l)  IR= 1 
GOTO  5 

17  IF( IR.LE.2) IR=2 
GOTO  5 


2j  print  ♦, 'Enter  both  strip  widths 

READ  *  , W  IN 

w(il=AMAXl(wIN(l) , W IN ( 2) ) 
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WT2)=AMINl(WlN(l),WlN(2>) 

!(?  =  ■» 
corn  s 

o 

c 

35  I=I-NE-1 

G0TP(2*3,3l,4-,5.,6.  i  o«  *  7  i  95  )  i  I 
G  IDGT  IS  AN  INTEGFP 

35  RRINT  ','tMtn  DK-tTS  OF  RESOLUTION  OtSIPEO 
READ  *,I0GT 
GOTO  5 

38  RRINT  , 'ENTER  N 

°£  A  0  *,N 
GALL  ISETN(N) 

I°  =  T 

so ro  5 

C  GO 

SC  3E  TUCN 

G 

C  PLOT  REOUE3T 
-.C  IF  (  IP  .NE..  I  GOTO  9  C 

IP  =  4 
R£TUPN 

S  PRINT  REQUEST 
63  IF(ie.NE..)GOTO  9C 

I  prR 
»t TURN 
50  STOP 

7,  I  F  C  I P  )  7  2  , 7  1 

71  RR IMT  9*3 

9*1  FOcMAT  ( *  Thl  IMPEDANCE  OF  TWO  Tn  IN  0  A  KA  LL  c.  l  STPIR?  C  t  N  Tt  RED  WITH  P 
llSPEGW*  TO  THEIR  MlLLINF  IS  CALCULATED.*/ 

,3*  RNTE-’  PAvAMlTEP  OP  FUNCTION  KEYWJRD  FROM  DElCW  .1ST  TC  0PECATE. 
*•/«  (U3UAi.LV  THE  FIRST  TWO  ChAP  v,F  DESCRIPTION)*/ 

„»  ,0  KEY  WU  R  0  INDICATES  an  OUTPUT  VAPIARlE*/ 

C*  RESOLUTION  IS  OPTIMIZED  WHEN  ICoT=;,  MAY  GE  'EDJGED  TD  I03T*/ 

C*  DIGFT;>  IF  IDGT  IS  NuN-Z ER 0 . . . . ♦ / 

C-  OISPlAY  wIlL  PRESENT  CURRENT  VALUES  */ 

S) 

I c  =  1 

’2  °  R I N  T  95  I , ( Q( I ) ,Kw ( I) ,UN( II  .  1=1 . L> 

95.  format  (  Of  S'"  -  IPT  I  UN  KEYWORD  ,  T  2 1  , '  JN I T  S  •  /  25  (  1  *  ,  A 1  .  ,  T 1 5  ,  AS  ,  T  22  ,  A  i  Z 
./)  I 
GOTO  5 

Mo  RRINT  951, UC  >  .FLOAT  (Nl  ,1H  ,  ( Q  (  I  ♦  1 )  ,  V  <  I)  ,  UN  (  I  +  l )  ,  I  -  1  ,  9  E  I  ,  G  (  9  i  ♦  3  )  , 
Cr  LOA  T ( IDGT  >  ,  1H 

951  FORMATf*  (7ARlMULfc',Tl2,"t/ALUE  ,/15(lX,AlS,19=,lPjll.(*,<.Y,AiC/)) 
GOTO  5 

93  ’PINT  *,*NC  OUTPUT  WITHOUT  A  GO  !!• 

GOTO  R 

C  SA  V-  Rt'uUrST 
95  IF ( IR.Nt. C ) GOTO  9  G 
w  P I TF ( ? , 96  . ) 

96-  TOPMAT ( 1H1 ) 

WPITT  ,95:iam,FLuAT(N),ln  ,(Q(Iy1),V<I),UN(I»1»,Ir1,NE»,0(NL+2) 
I, FLOAT (IDGT) , 1 h 

IP=f> 
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RETURN 

ENO 


subroutine  curr<n, isa,suma,sumb, iso 

CURRENT  IN  ELEHENTRY  STRIPS  FOR  THE  HALF  PLANE  SYMMETRY  CASE 
ARE  CALCULATED.  ARRAY  ALLOCATIONS  ARE  MADE  HERE  FOR  3lAN<  COMMON 
IN  THE  MANNER  DESCRIBED  IN  FUNCTION  INA,  OTHER  SUBROUTINES 
REQUIRING  SPACE  MUST  OBSERVE  THE  NOTES  BELOW. 

" r«~rS~A“¥nT5*~SP'£ClT INC  THE  MATRIX-  NUHTSTR  ^FOR  VA  & IA3LES  AS  FOLLOWS! 


m  mjtual  inductance  hatrix  is 
phi  matrix 

PSI ( INVERTED  I  MATRIX 
PSI  MATRIX 

MOW?  ARTS  'MATRIX - 


INA(2,ISA(1) 
INA  (2  ,  IS  A  ( 2 ) 
I NA  <  2  *  IS  A ( 3  ) 
INA  (2«ISAI4) 


THE  WORK  AREA  EXTENDS  6»N  BEYOND  INA (2 ,ISA (5 » ,N2 ,  N2 I  AND  IS  PROVIDED 
BY  3  2*N  VECTORS  OEFINEO  BELOW,  IF  N  VECTORS  ARE  DEFINED 

ELSEWHERE,  THEY  WILL  BE  HP  ITT  EN  OVER,  THE  AREA  MAY  BE  USED 
FOR  VOLATILE  STORAGE. 


IT 
C‘ 

c 
c 
c 
c 

~C - 

C  TWO 

e  the 
c 

*  VECTOR 
B  VECTOR 
CvECTO'R- 


INfl(2,ISA(b) 


2*N  WORK  AREA  VECTORS 
THlRO  2*N  AREA  IS  USED 


C 

C 


FOR  1/R  INV( 2  , ISA ( S ) 

INV(2  ,ISA(TI 
I NV  (  2  ,  IS  A  (  8 ! 


TNVt  2 ,TSA(9J 


EXTERNAL  INA.SZINO 

DOUBLE  PRECISION  SPSI , SPS IPH , SPS IR, SPS IPMR, DENOM 


COMMON  //A ( 1 ) 


COMMON  /STR IPS/FREQ, SEP, W( 2) ,EW(2) ,T  HK ,SKD°S, £TH< 
“C7TrHO,EOC,Paw,IOGT,  IER,PI2 


C 

C 


DIMENSION 

DIMENSION 


RESOC ( 2 ) , NOI A (4 1 
I5AT9T  - - - 


C 

C 

C 

C 

t r 
c 

c 


THIS  VECTOR  DESCRIBES  THE  ATTRIBUTES  OF  MATRICIE  S  FOR  THF  IMSL 
OIAGONAL  MATRIX  TIMES  A  FULL  MATRIX  FUNCTION. 

OATA  NOIA/4*0/ 


BOO 

C 

c 
c 


PRINT  «,'*CURR* 

SET  MATRIX  LENGTH 

N0IA<1»  =NOIA<2)=N2=2-N 
FILL  THE  MUTUAL  INOUCTANCE  MATRIX 
TREF=SECONO(CPI 

.TTGOTO  T 

CALL  HUT (A(INA(2,ISA(1 1 ,1,1) I , A ( INA ( 2 , IS  A ( 21  ,1,111  ,N2,2*N2,SZlND) 
CONTINUE 

TMUT=SECONO<CP>-TREF 
IFIISK.EQ. 01  PRINT  9Cu,6HMUTUAl,TMUT 
FORMAT (IX, A8,*  SOLUTION  TIME=  *,Ffa.2,*  SECONDS*/) 


COMPUTE  PHI=H*M/R 

STUFF  A  BAND  ARRAY  WITH  1/R  USING  THE  SMALLER  OF  THE  TAPE 


cIoooom  i  oejooo  ,  noooq  ooo 


C  THICKNESS  OR  A  MULTIPLE  OF  SKIN  DEPTH. 

HF  =  PI 2*F REQ 

C  HOJlFr  US  ( SEL ATI VE  PERMEA BIL I T Y )  IF  ME  CONDUCTOR  IS  “ACNETIC  MATERIA 
UR=  1 . 

SK0P=SQRT<KH0/ <FR£Q*1 . £-. 7*UR> ) /P IE 
ETHK=“AHTN1  (THK, SKUPS *SKDP) 

IF(ETHK.NE. THKJPRINT  ♦, 'EFFECTIVE  TAPE  THICKNESS  AOJUSTED 
00  13  1=1,2 

lu  SESOC ( I) =RH0/(£W< II ♦ ETHK) 

c 

L0C=INV(2*ISA(6) ,1,1)-1 

LOD?=LOC*N 

RINV=1./RES0C<1) 

RINV2=1./RESDC<2) 

00  12  1=1, N 
A<L0CH)  =RINV 
12  A  <L0C2+I)=RINV2 

C  ^RTNI  r ♦TND  INr  - 

C  SALL  MATP(A (INA(2,ISAI1I ,1,1)1 ,N2» 

c 

C  PRINT  * , *  *  1 /R  ♦, (A  (L0C*I>  ,I  =  1,N2) 

c 

C  COMPUTF  PHI  DO  H*M  FIRST 

- C ALT  SHUL“H CAT TNA  C2,  ISA  f  1)  ,1,1)1  ,  fi  t  INA  t?,rSAT4M  IT, ITT  ,N2  ,WFT 


PHI  =W*M/R 

CALL  VMULFB(AIINA(2,ISA(4) ,1,1)) »N2,A(L0C*1),N2,NDIA,A(INA<2,ISA<2 
C) ,1,1)1 ,N2) 

. — PRTNT  '  M  r*PHr  a"  " 

CALL  MATP(A(INA(2,ISA<2) ,1,11) »N2 ) 

COMPUTE  PSI( INVERTEO) 

NOH  00  CW*M/R)*W*M 

CALL  VMULFF  (A  (INA  (2, ISA  (2)  ,1,1)  )  ,  A  ( I  N  A  (2  ,  I SA  (  *  )  ,1,1)  ),N2,N2,N2,N2  • 

- — (CNZTBTTTtA-rrrrsA  far, i,D)  ,n2,ier) 

IF(IER.NE.O) PRINT  * ,  'ERROR  IN  PSIUNV)  MULTIPLICATION,  IER=  ', 

CMOO  C IER, 32 ) 


PRINT  PH  I ' W*  M  • 

- TTAETTNIATPrATTNAC?,  ISA(3)  ,1  ,trr  ,NET  -  -  -  - 

ADD  R 

L0C=INA(2,ISA(3I ,1,1)-1 
L0C2  =  INA  (2, ISA  (3)  ,N*1,NMI  *1 
00  2C  1=1, N 
J=I*< I-1)*N2 

— — ffrcTjc>  jr=A  {Loc+jr^REsoca)  - - 

3  A<L0C2* J> =A<L0C2*J) +PES0C < 2) 

PRINT  *, **PSI(INVERTED>=  * 

CALL  MATP(A<L0C*1I ,N2) 

INVERT  THIS  TO  OBTAIN  PSI,  THIS  DESTROYS  THE  1/R  VECTOR 
TTS  TNE~  EXTENDED  WORK  AREA  IS  USED. 

I0GT2  =  I DGT 

CALL  L IN V2F { A ( INA ( 2 , ISA (3 ) , 1 , 1 ) ) ,N2 , N2 , A (IN A ( 2 , 13 A ( 4 ) , 1 , 1 ) ) , I0GT2, 
CA< INA  <2, ISA (5) ,1 ,1) ) ,IER) 
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IHIE^tNE.  J ) PRINT  •  ‘ERROR  IN  INVERTING  PSI(IMV),  IER=  ', 
C*tOO(  IER, 32) 

PRINT  *  *  *  *  I OGT  IN  INVERTING  PSI*  ',IDGT2 
°RINT  »,'*PSI  * 

CALL  “A TP (A ( INA <2, ISA <41 ,1,1)1 ,N2) 


NOW  00  PSI’PHI  INTO  WORK  AREA 

CALL  VMULFF<A< INA  <2, ISA (41 , 1 , 111  , A < I N A ( 2 , 1 SA C 2 > , 1 , 1) ) , N2 ,N2, N2 , N2 , 
CN2, A (IMA <2, ISA <5 1 ,1,11) ,N2,IER) 

IF(IER.NE. Cl  PRINT  *  ,  'ERROR  IN  PSI’PHI,  IER*  ' , MOOC IER, 32) 

PRINT  *,'*PSI*PHI=  * 

_ '  :~tll  hutpt a  sttstti;  nr  ,N7i - 


PERFORM  SUMS  OF  MATRIX  ELEMENTS  FOR  E  AND  F  SOLUTIONS 
SUM  »SI 

CALL  SUMAL'LTirriNAI?,  ISATU^r*  ITT  , N2,5FSIT - 

SUM  PS  I*  3  H I 

CALL  SUMALL(A(INA(2,ISA(5) ,1,111 ,N2, SPSIPH) 

SUM  THE  RIGHT  SIOE  OF  PSI 

CALL  SUMRT (A  < INA (2, ISA  (4)  , 1 , 1 ) > , N2, SP SIR) 

SUM  RIGHT  SIOE  OF  OF  PSI*PHI 

CALL  SUMRT  (A  ( IFTAT2 , 1  SA  <51 , 1 ,  IT  )  ,N?,SPSTPHR1 - 


DENOM=SPSI*’2+SPSIPH**2 

E*-(SPSI’(-SDSIR)*SPSIPH* (-SPSIPHRII  /DENOM 
F=(SPSIR-E*SPSI) /SPSIPH 

PRINT  ,  'V0LT1TGT  DROP  IN  BOTTOM  (WI3E)  TAPE*  ‘r,TP: - 

PRINT  *,  '*E1=  * ,2. ’SPSI 
PRINT  '  ,  *'  £2=  »,2.*SPSIR 
PRINT  *,'*£3=  ',2. ’SPSIPH 
PRINT  , '*E4=  *  ,  2,  ’SPSIPHR 

"  TTTAT*  S 1 C  OWDTC’PT^TREF^TMUT - - 

PRINT  9JQ,7HV0LTAGE,TEAF 

SET  UP  SOLUTIONS  FOR  A  AND  B  IN  THE  FORM  OF  IMSL(AX=B>  OR 
PSI(INVI*A*<E*PHI’F> 

PUT  RIGHT  SIOE  IN  SOLUTION  VECTOR 

CALL  V M U'L M <  ATT  NAT2 ,TSTT 2T,1 , IJ  ) ,'F , N2VffTTWr? , ISA r7T,T,~TrFT 
L0C=INV<2, ISA <71 ,1,11-1 
L0C2=L0C*N 
JO  3  j  I  =  1 » N 
A  (LOCH)  =A  (LOC  +  I)  ♦£ 

A (L  0C2 ♦ I )=A ( L0C2+ I ) *E-1, 


JO  SOLUTION 

I OG  T2  =  TOGT 

CALL  LEQT2F < A  < I N A  < 2 , ISA <3  I , 1 , 1 ) I , 1, N2 , N2 , A <L 3C ♦ 1 ) , I0GT2 , A  TINA (2, IS 
CA  (5  )  ,1,1)1  ,  I  ER) 

C  “  “PRINT  A“  STJL"  =  '  ,  (  ATLOC+  I )  ,T=X,TJ7T - 

IF< ICR.NE. Cl  PRINT  *, 'ERROR  IN  A  SOLUTION,  IER*  • , MOD < IER , 32) 

PRINT  *,'»IOGT  IN  A  SOLUTION*  ',10312 

C 

38 


s 


C  NOW  QO  FOR  8 
C  PSI<lNV)*8=<PHl*E-F> 

C  CALL  MATPJA (INA (2, iSA<2) ,1 ,1) » ,N2) 

CALL  VMULMP< A  ( INA  (  2  ,  ISA  ( 2 )  ,l,l)),E,*2,A<m<2,ISA<8),l,im 
LOCB=IN\f (2* ISA (8 > 

DTT4j_T=T7N2 - - - - 

40  A(LOCB*I»=A(LOCB*I)-F 

I0GT2=IDGT 

CALL  LEQT2F(A(INA(2,ISA(3) , 1* 1)1  , 1,N2 ,N2 , A < LOC B*l>  , 193 T2 , A ( I NA ( 2, 
CS A ( 5) *  1 » i) > ,  IER) 

C  PRINT  ♦,**a  SOL=  *,<A<LOCB»I),I=l,N2> 

IFTIFN.  NE.O)  PRINT  *T*ERROR  IN  B  SOLUTION,  IER  =  *,M0D(IER,32) 

PRINT  ,  •* I OGT  IN  B  S0LUT10N=  *,10312 
TSOL 1= SECOND t CP) -TREF-TEA F-TMUT 
PRINT  9oO,/HCURRENT  ,TSOLI 
C  SUM  REAL  ANO  QUADRATURE  CURRENTS 
C  A  =  REAL,  9=  QUADRATURE  CURRENT,  C=  MAGNITUDE 
Locc=iN/rzyrsjTFjTTr,  i  y-i 

LOCC2=LOCC«-N 

SUMA=S'JMB=0. 

SUM A2=SUM82=0 
DO  St  1  =  1, N 
SUMA=SUMA«-AIlOC*I) 

atlocoii = a  (L'oc*it*a  (Locrrrr+ a  CLOCB>n  *ailocbhi 

AluOCC2*II=A<LOC4l*N)*A<LOC*I*N> ♦  A <l OCB* I* N) *  A  (LOCB*I»N) 
SUMA2=SUMA2*AIL0C*I»N) 

SUMB2=SUM82» A  (LOC8«-I*N) 

80  SUMB=SUMB*A(LOCB*I) 

SMAG1  =  SQRT  tSUMA**  2  +  SUMB**2 ) 

- SNAGZ=SQRTTSUMA2*’2*'SUMB2**2J 

SUMA=2.*SUMA 
SUMB=2.*SUMB 
00  52  1  =  1, N 

A<LOCC>I)=SQRT ( A  <  L  OCC ♦  I) ) 

52  ATLOCC*I*NJ  =SQRT IA (LOCC  +  I *N) ) 

~P  R  TNT  *  ,  IFTER  E  NC  E  IN  CURRENTS  =  • ,  A9SISNAG1  -SMA32)  f\ 3M  AG1  *SMAG2 

C* 2 OC .  ,  *  PERCENT* 

C  PRINT  ♦ , ** SUMA= • ,SUMA, •  SUMB=*,SUMB 

RETURN 
END 
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SUBROUTINE  VMULM(A,V,N,0) 

C  SPECIAL  PURPOSE  SQUARE  MATRIX  TIMES  A  VECTOR  WHOSE  ELEMENTS 
0  ARt  IDENTICAL 

DOUBLE  PRECISION  Sl,S2 
DIMENSION  A <N,N> ,0(N) 

00  2GC  1=1, N 
S1  =  G. 

00  100  J=1,N 
10d  S1«S1+A(I«J> 

200  0(II=S1*V 

RETURN 

CNfRV  "VMULMP 

C  SUBTRACT  one  from  the  last  N/2  vector  elements  AND  MJ.TIPir  as  above 
NQ2=N/2 
00  160  1  =  1, N 
S1=S2=0. 

00  ISO  J= 1 , N02 
S1=SI+A ( I, J) 

15  li  S2  =  S2«-A(I,J+N02) 

160  0(1)  =V*  (Sl«-S2» -S2 

RETURN 
ENO 
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SUBROUTINE  SMULM < A  ,B,N.S> 

C  SPECIAL  PURPOSE  SQUARE  MATRIX  TIMES  A  SCALER 
OIMENS ION  A (N,NI «B(N«N) 

00  100  J=1*N 
00  IOC  I =1 » N 

Tiro - «m,ji=A  lit  jj*s 

RETURN 

C 

ENO 
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SUBROUTINE  SUMALL ( A »N, S ) 

DOUBLE  PRECISION  S,SUH 
DIMENSION  A (N»N) 

C  SUM  ALL  ELEMENTS  OF  A  SQUARE  MATRIX 


C  SUM  ALL  ELEMENTS  IN  THE  RIGHT  SIOE  OF  SQUARE  MATRIX 
ENTRY  SUMRT 
JS=N/2*1 
SOTO  199 


ovnoooovfi  I  -r  o  wo  |  oo  o  o  o  o  o:  o  o  o  o 


FUNCTION  SZIND(PP,EP,EwP,IER> 

THIS  FUNCTION  RETURN  THE  MUTUAL  INOUCTANCE  FOR  THIN  STRIPS  OF  WIDTH 
E  HP  METERS  SEPARATED  BY  PP,  ANO  DISPLACED  BY  EP  USING  THE  METHOD  OF 
BROOKE  ..  ET  AL  IN  'INDUCTANCE  ANO  CHARACTERISTIC  I MPE D ANC E . . . ' 

25  OCTOBER  1966,  (APPENDIX)  WITH  THE  LENGTH  TERM  REMOVED. 

OUTPUT  IS  IN  UNITS  OF  MICRO-HENERYS  PER  METER*  THEN  *1.E-C6 
LIMITS  ON  CALCULATION  ARE  PROVIOED  WITH  AN  ERROR  FLAS,  I ER  *  TO 
INSURE  THE  INTEGRITY  OF  THE  ANSWER,  CONSISTENT  WITH  THE  COC 
660:  COMPUTER.  (29  SIGNIFICANT  DIGETS  IN  THE  CALCULATION) 

DOUBLE  PRECISION  A 1 ,81  ,C  , D  , P , E  ,E  W .DEL. 

DIMENSION  EHP ( 2) 

DATA  STD/1. Eil/ 


TST=SQRT(PP*  2+EP* -2) /ABS ( EWP) 

I F ( TST .LE. STD) GOT 0  4 

IFdER.EU.  0  )PRINT  951 , EWP",~PP,  EP ,  TSTYS TO 
51  FORMA  T ( *  ERROR  IN  MUTUAL  INOUCTANCE  C ALCUL A T I  ON* / T 1 0 , 

C  A  WIDTH  OF*  »  G 1 0  •  4  »*  METERS  AT  SPACING  OF*  ,  2 ( 1 X , C 1 0 . 4, H , > , 
C/T10,*YIELOS  A  SPACE  TO  WIOTH  RATIO  OF*,S10.4,*  EXCEEDING  THE  *, 
C* ALLOWABLE  OF*, GIG. 4/) 

TER  =  i E  R  + 1 
RETURN 
CONTINUE 
P  =  PP 
E  =  EP 

EW=EWP(1) 

OELSITWP  < 1) -EWP ( 2  ) 

A1=B1=C=D=0.0 

IF(PP.EQ.O..AND.EP.EQ.O.) GOTO  91 3 
IF ( PP.EQ.  u . 0)  GO  TO  90 G 

A1  =  P  <(EW-E)  DAT  AN2  (  (E-EW) , P) ♦ < E- DEL >  *DA TAN2  (  C E-0 EL )  , => )  * 

CE’OAT AN2  (E , P ) - ( E* EWP (2 ) ) * OA T AN2 ( <£*EWP(2>)  ,  =  )  ) 

TO  TO  901 
90  C  CONTINUE 

IF ( E.EO . Ew ) GOTO  902 

901  31= (P*P-(E-EW) **2)*0LQ&(P*P*(£-EW>**2)/4. 

902  CONTINUE 

C=(E*E-P*P) *DLOG (P*P*E*£> /4. 

C  ♦" f  ( E  -  O  EL  )  M  E-  O  EL  )  -P*~Pr*DrOG  (P*P*T  EOED  «  TE-OEL  >  )  /A. 

903  CONTINUE 

0=(P*P-(E*EWP(2) )**2)*DL0G(P*P+(E*EWP (2) )**2) /4. 
ANS=(A1*B1+C*0)*0.2/(EW*EWP(2)  ) 

GOTO  53 

10  A  NS  =  0 .2*OLOG(l./EW) 

PRINT  * , **SELF ' 

CHANCE  TO  UNITS  OF  HENERYS 

0  SZIN0=ANS*1 .E-06 

PRINT  *, 'P,E,EW,ANS  *,P,E,EW,ANS 
PRINT  '  ,  * A1 ,B1 ,C »D  ‘,A1,B1,C,0 
RETURN 
END 
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SUBROUTINE  MOX ( I  * J , II, J J» 

C  CONVERTS  A  REFERENCE  <I,J)  TO  ELEMENT  NUMBERED  PER  BRDOKE  ET  A., 

C  TO  EQUI VALANCE  ASSUMING  HALF  PLANE  SYMMETRY.  CONSEQUENTLY, 

C  2*<N**2+2J  MUTUAL  CALCULATIONS  ARE  REQUIRED,  RATHER  THAN  16*N**2 

C 

~c - - - — 

COMMON  /STOR/N,MAXN,ICCRE,NV <6»,NA<6) ,MM(6) 

C 

C  THIS  FUNCTION  MILL  TEST  0  IF  LOWER  TAPE,  1  IF  UP3ER 
ITOP  OO  = <MOO(K-i ,2*N> >/N 

C 

X”XtJT  SWALtER  IN  II - 

1 1  =  MIN J ( I, J» 

JJ=MAXQ (I, J) 

C  LOCATION  FLAG  FOR  I 
IT  =  ITOP (II) 

JT  =  I T  OP ( J J) 

X  TH IS~FL  A G'^W ILL  TEST  *  IF  J  J  TOP  ANO  If  "B0TT01,  -  IF  "TT'TDPy  JJTTOTTOH  . 
I0P= JT- IT 

C 

C  SET  THIS  MULTIPLIER  TO  2  IF  BOTH  ARE  TOP,  1  OTHERWISE. 

IH=IT*JT*1 

C 

rFTTOPT5“,  113,10 
C  I  TOP,  J  BOTTOM 
5  JTEHrJJ 

JJ=II*2*N' 

II=JTEM-2*N 
RETURN 


L  Itil  LE  li.  IN  Ltr  l  MALE  PLANE 

1 G  IF( II-IM*N) 15,15,17 

C  LEFT 

15  IF(IOP»21,21,16 
C  LEFT,  II  BOT,  JJ  TOP 

16  RETURN 

C  BOTH  IN  RIGHT  HALF,  PUT  IN  LEFT 

17  II=II-?*N 
JJ=JJ-2*N 
SOTO  15 

C  BOTH  ARE  TOP  OR  BOTTOH 
21  IO=II-IT*N-l 

- IT  =  II_I0 - 

IF< JJ.GT.2*NIGOTO  25 

JJ=JJ-IO 

RETURN 

25  JJ=JJ*IO 

I F  f JJ.LE.N* (2»IM> >  RETURN 

RETURN 

END 


SUBROUTINE  NUT ( AN , A W,N2 ,N4 , FUN > 


STUFF  A  2N*2N  MATRIX  WITH  MUTUAL  IMPEDANCE  USING  FUN» 

REPRESENTING  A  4N*4N  MATRIX  WITH  HALF  PLANE  SYMMETRY,  THE 
LONER  HALF  IS  IGNOREO  AND  THE  UPPER  HALF  FOLDED. 

EXTERNAL  FUN 

FUN  IS  THE  IMPEDANCE  CALCULATION,  £w(2>  CONTAINS  THE  INCREMENTAL 
STRIP  WIDTHS  FOR  LOWER  (THE  WIOER)  AND  JPPER  TAPE  RES3EC  T IrfELT . 

DIMENSION  AM(N2,N2) ,AW(N4,1) 

DIMENSION  EWB(2) ,EWT<2) 

AW  IS  A  WORKING  SPACE,  MUST  BE  (2N)**2*4N  LONG,  MAT  BE  SJPPlIED 
BY  A  2N*2N  ARRAY,  BUT  CAUTION  THAT  THE  NEXT  4N  LOCATIONS  WILL 
BE  USEO. 

MM,*****************************-***#**************** 

NOTE  THAT  STORAGE  IN  AW  IS  ROW  AND  COLUMN  REVERSED 

COMMON  /STRIPS /FREQ, SEP, W(2I , EW ( 2) , T HK , SKDPS,  E TH< 
C,RHO,EOC,POW,IQGT,IER,PI2 
PRINT  * , **MUT  * 

RESET  ERROR  COUNTER  IN~FWcTTO N“FUN  ("MUTUAL  INDICT A  N  C  £  SUBROUTINE) 
IEP  =  3 

COMPUTE  ELEMENTARY  STRIP  WIOTHS 
OD  1  1=1,2 
EW(IJ=W(I)/N2 

ToaqT"  widths  pAIrs~for  fun"  calls 

DO  2  1=1,2 
EWB  (II =EW( U 
EWT(II=EW(2) 


N=N2/2 

N3=N2*N 

OFFSET  OF  TOP  SHORT  TAPE  FROM  ORIGIN 
El=(W(ll«W(2>>/2 

OFFSET  OF  RIGHT  EDGE  OF  TOP  TAPE  FROM  ORIGIN 
E2=E1*W(2) 

STUFF  THE  AW  MATRIX 

=>RINT  *,  **AW  ELEMENTS  *  ,4 *  ( N** 2«-N  ) 

DO  1,  J= 1  ,N 

AW(J,1)=FUN(3,EW(1I*(J-U ,EW8, IER) 

AW( J4N,1»=FUN(SEP,E1*EM(2)*( J-l) ,EW, IER) 
AW(J+N2*l)=FUN(u*N(ll-EW(ll*(J)*EW3*IER) 
AW( J»N3,1> =FUN(SEP,£2-EW(2)*( J) ,EW,IER> 
AW( JtN,N»l> =FUN(0,EW(2I* ( J-l) ,EWT,IER ) 
JW(  J»N3,N*1)=FUN(4,EW(2)»  (N2-JI  ,FWT,  IER) 
PRINT  ", *N#S=»,N,N2,N3,N4,NPl 
DO  30  1=2, N 
EO=EW (1) * ( I - 1 » 

00  33  J=  1  ,N 


,  '  ' 


OOOl-  ooooo 


4WC J*N, I)=FUN<S£P,A8SCE1*EW(2I*< J-l> -E0> .EW.IER) 

30  4W< J»N3,I)  =FUNCSEP»  E2-ENC  2IMJ)-£0,EW,IER) 

NP1=N*1 

PRINT  *,**AW  PRINT  LIMITS*  * ,N4 , NPl , < (AM < J, I J , J=I , N4) , 1= 1 , RP1 ) 
PRINT  '  ,  *  A  M  ' 

NOW  STUFF  AH 

30  1C  J  11  =  1. N2 
00  ICO  JJ=1,N2 
CALL  M0X(II,JJ,I,J) 

CALL  M0X(II,N2»JJ,IA,JA) 

00  AM(II,JJ)=AH(J,II ♦AWCJA.IA) 

PRINT  *,  (CAM (I, J»  « J=1,N2»  .I  =  1,N2) 

PRINT  * . *  *MUT  DONE* 

1F(IER.EQ.0)RETURN 
PRINT  950,  IER 

950  FORMAT ( *  MUTUaL  INDUCTANCE  ERROR  COJNT=  *,I7> 

RETURN 

END 
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^UNCTION  INA(M,K,I,J) 

C 

C  THIS  FUNCTION  PERMITS  SPECIFICATION  OF  ANY  (I,J)  ELEMENT  IN  ANY  (K) 

C  SPECIFIED  (N*M)**2  ARRAY  OR  N*H  VECTOR#  WHERE  M  MAY  3E  1,2,4,8,12,16. 

C 

C  ST ORACTE  13  In  3 LANK  COMMON,  ARRAY- 5TIC0RE)  NTT-PTHE  SQUARE  ARRAYS  F  I 
C  FIRST.  THEN  THE  VECTORS. 

C 

C  THE  NUMBER  OF  EACh  ARRAY  TYPE  QEFINEO  IS  STORED  IN  NA(IMI  AND  NV(IM) 

C  '-MERE  IM=1,6  (SPECIFIC  THE  M  IN  ORDER  SIVEN  ABOVE) 

C 

COMMON  //All) 

COMMON  /STOR/N,MAXN,ICORE,NVI6>,NA(8> ,MM(S) 

DIMENSION  I  £  NT  I  2 ) 

C 

DATA  IE.NT/bHVECTOR,5HARRA  Y/ 


C  ENTRY  FOR  ARRAY  INDEXING 
IT=2 
JJ= J-l 

C  SET  INDEX  OFFSET  TO  ARRAY  NUMBERS 
NO=  6 

ix  =  : 

c 

c 

c  SEE  IF  M  VALID 
3  DO  5  IM=  1 , 6 

IF(M.£G.MM(IM) )GOTO  1C 
5  CONTINUE 

7  PRINT  *,'CAlL  to  inoex  error,  CA  LL-  *,M,«,I,J 

print  * , IENT (IT)  *  *  ENTRY  POINT* 

I NA  =  1 
RETURN 

C  INDEX  TO  NO  OF  M  ARRAYS  OR  VECTORS 
13  CONTINUE 

C  GET  ADDRESS  OF  FIkST  M  ARRAY  -1 
IF(IM-l) 12, 3C 
12  L=IM-1 

00  IS  IM=1,L 

15  IXrlXYNV (NO+IM) *(N*MM(IM) ) **IT 

C 

C  CHECK  FOR  VALIO  K 

SC  IF(K.LL.O.OR.K.GT.NV<No+IM) JGOTO  7 

c 

C  CHECK  FOR  VALIO  J 

IF( JJ.LT.O.OR. JJ.GT.N*M)SGT0  7 

c 

C  CHECK  FOR  VALIO  I 

TF(I.LE.O.OR.I.GT.N*M)GOTO  7 

C 

C  SET  REQUEST  INOEX 

INA  =  IX»(K-D*  (N*M)  -'.T  +  JJMN’MMI 

RETURN 
C 

C  ENTRY  FOR  VECTOR  INOEX 
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cntry  inv 
I  T=  1 
J J  =  NQ  =  0 
I X= IN2 
SOTO  3 
T 
C 

ENTRY  ISETN 
M=MAXC<M,2I 

C  THIS  ENTRY  COMPUTES  A  MAX  N  BASED  UPON  MA,  NY  AND  COMPARES 
C  THE  N  REQUESTED  TO  IT,  SETS  TO  MAX  IF  MORE  REQUESTED. 

C“ 

C  SUM  THE  NUMBER  OF  N**2  SUBARRAYS 
NN2A  =  0 

DO  10-  I M  =  1 » 6 

103  NN2A=NN2A»NA( IMJ *MM( IMI **2 

C 

IT  SUHTHE  NUMBER  OF  N  VECTORS 
NN  V=  0 

DO  110  IM  =  1 ,6 

lie  NNV=NV< IM» *MM( IM) ♦NNV 

C 

C  QUADRATIC  SOLUTION  FOR  MAXN 

- HTVNrrSQRT TFTNV* k .’•rNR28*TCORE7  -NWT  /T2.*NN2A  ) 

IF(NN2A.EQ.G>MAXN=ICORE/NNV 
115  TN2=NN2A*N**2 

IUSE0=NN2A*N**2»NNV*N 

PRINT  *,'N  MAX=  '  ,  MA  XN  »  '  UNUSEO  ODRE*  I  CO RE-I US  ED 

IFtM.ST.MAXNIGOTO  200 

- R=H - - - - - —  ■ 

INA  =  N 

C  NORMAL  RETURN 
RETURN 
C 

C  ERROR 

Tnj - ICN=NNV*N*NN2A«N*N - 

N=MAXN 

INB=N 

PRINT  *,•  ERROR,  REQUEST  FOR  SUBDIVISION  TDD  .ARSE,  SET  TD  *,MAXN 
C,'  CORE  REQUIRED  =  * , ICN 
SOTO  115 


3U0ROUTINE  MATP(A,N2> 

OIMENSION  A  »N2  ,N2  » 

DIMENSION  I F  OR ( 2 ) 

DATA  IFOR/1H  ,luhlX,C>10.4)  J  / 

C  LIMIT  PRINT  WIDTH  TO  CARRIAGE  SIZE 
NC=mIN3 (i3,N2) 

ENCODE <10,000, IFOR) NC 
Djj  FDRMA T <1H< , 18, lh ( ) 

°RINT  IFOR, ((ATI, J),J=1.M2»  *  I=1,N2) 

RETURN 

END 
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SUBROUTINE  PLOTI < A,H,N2,W) 

C  4  CONTAINS  30TT0M  TAPE  CURRENT  (MAGNITUDE)  IN  THE  FIRST  N  LDCS, 

C  AND  THE  TOP  TAPE  IN  THE  lAST  N,  THE  TOP  TAPE  IS  SMALLER  OR 
C  EQUAL  TO  THE  TOP  IN  WiuTH.  OUTBOARD  ELEMENT  FIRST . 

C  H  CONTAINS  THE  TAPE  WIDTHS,  LARGER  FIRST 

c 

C  AN  ALLOWANCE  FuR  H  SHOULD  BE  MAOE, 

C  THE  FOLLOWING  METHOO  USES,  AS  A  MAX, FOUR  N2**2  AND  SIX  N  ARRAYS 

C 

c 

DIMENSION  A(N2),P(1),SC(2)  ,  W  ( 2 ) 

EQUIVALENCE  (AMX,SC(1)» 

C  THc  FOLLOWING  IS  FOR  Th£  PLOTTING  ROUTINES 

c 

COMMON  /CKT3UF/IPP ,M3UFF( 1C2A1 
DATA  IO3/1024/ 

DATA  SC ( 2) / G . / 

DATA  IP/3/ 

C 

IF  ( 1° ) 2 , 1 

1  IP=1 

'  "  HR  I  NT  RlO 

HID  FORMAT{*  PROGRAM  PAUSES  WITH  CURSOR  AFTER  EACH  PLOT,  ENTER  A  CHAR 
CTO  CONTINUE*/*  AFTER  THE  SECOND  PLOT,  ENTRY  DF  THE  FIRST  LETTER  OF 
C  THE  3£  LOW  DESCRIPTION*/*  WILL  CAUSE  THAT  RESULTI*/ 

CT5,*REPLOT  BOTH  TAPES* /T5 ,* FULL  WIDTH  SMAL.ER  TAPE  *  REPLOT* / 
CT5,*PL0T  SIZE  CHANGE  (TWO  CURSOR  INPUTS  DEFINE  FRAME  SIZE)*/ 

.  ■■C1-  -  - - - - - - - - - 

2  N=N2/2 

ISPCr**  "N2*N2«-3*N2 
ISPC=MINu( ISPC,199I 
C  FIND  THE  LARGEST  VALUE  FOR  SCALING 
A  M  X  =  0  . 

TO'  It  T=IYN  - - - - - 

1G  AMX=AMAX1(AMX,ABS(A(I)I) 

C  SCALE  VERTICAL  SCALE  TO  CURRENT  DENSITY 
NPl=N*i 

DO  12  1=  NPl  ,N2 

12  AMX  =  AMAX1(  AMX,ABS(A  III*X(1I/X(2HI_  _  _ _ _ 

C  Plot  LOWER  TAPE 
IS  DO  2C  1  =  1, N 

23  =>(I)=P(N2*1-I)=A(I) 

C 

C  SET  SCAIF 

CALL  PLOTXY(SC,2, IR,4) 

c 

C 

call  PLOTXY (P,N2, IR, 1) 

C 

c 

C  ATTEMPT  TO  CENTER  SMALL  TAPE  FOl?  MORE  PLEffSTNC  D I SPLAT-  AND  EXPARTI - 

C  TO  1 G 3  TO  2x0  POINTS 

A00=(N2-1)‘ (W(l)/W(2»-1. ) 

M=1 
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IF!A00.NE. u,IM=<ISPC-l>/!N2-l*ADD> 

M=MAXC!1,M» 

C  CHECK  TO  SEE  IF  UPPER  TAPE  SHOULO  BE  ZERO  FILLED 
IF  !ADQ*N2*M,GT  .ISPC1GOTO  22 
IF ( W ( 1 I /  H  (  2 ) .LE.16IG0T0  30 
22  CONTINUE 
AOO  =  3 
M  =  1 

IF  I W ( 1 ) / H ( 2 l .NE.1.)PRINT  902 

902  FORMAT!*  TOP  ANO  BOTTOM  STRIPS  PLOTTED  TO  DIFFERENT  wIOH 

C> 

C  ' 

0  STUFF  ZEROS  IF  REQUIRED 
30  NA=(A09/2.» ** 

NT  =  2  NA*  M* ( N2- 11*1 
DO  32  1=1, NA 
32  R<II=P! NT ♦1-11=0, 

C  MOVE  ANO  PLOT  UPPER 
35  NM1  =  N - 1 

DO  1=1, NM1 

oo  a:  J= 1 , M 

•*0  R(NAMI-1)*M*J»=P(NT  +  1-NA-(I-1)*M-J)  = 

,,!  <  1  A!NM>1)-A(N+I»  »*FLOAT  !  J-li  /FLOAT  (M»tA!N*in*A(l»/4(2n 
M2=M*1 

C  PRINT  *,**P  IN  * ,  IPII)  ,1=1 ,NT> 

C 

C  FILL  in  center 

DO  a2  1=1, M2 

♦  2  P(NAMN-1»*M*I*=A<N2»*W<1»/W!2) 

C  PRINT  -,*P  Ou T=  *,(P(I),I=1,NT) 

C 

n  3  CONTINUE 

CALL  PLOTXY  IP,NT,IR,1) 

IFUR.fq.xhRIGOTO  15 
IF (IR.EQ.IHPHjOTO  5 C 
IF( IR.tQ.lHFlGOTO  22 
C 

RETURN 

51  CALL  PlOTXY!P,N2,IR,2> 

RETURN 

END 


-UNCTION  TIMEL(TL» 

CALL  CPUriMe  (1ST) 

C  RRINT  1ST  .  FLOA  T (ANO( .NO T .MASK ( 35  »  ,SHIFT<IST.-36>)) ,  FLOA  T ( ANO( 

C  C NOT . M  A  3K ( 4b  I  ,  SHIFT  (1ST* -12)  n,O.Ojl*cLOAT(ANH.NOr.MASK<48*,ISTM 

FORMAT  ( IX  ,02u/*  Av/AILABLE,  SEC  NOdt  MILE  NOrf* /IK , 3 F6. 3 > 

TIMEl  =  TL  =  FLOAT  <AN0(.N0T.MASK<36>  .SHIFT  (1ST,  -36)  )  >  -  ( FED  A  T  ( A  >1  □  (  .  NOT 
CMASKl4b)  .SHIFT  (1ST  «  -  12  •  )  )  ♦O.uol’fLOAT  (ANO(.N3T  .M4S<(A9)  .1ST)  )  > 
RETURN 
E  NO 
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SUBROUTINE  PLOTXY ( Y,NP,IR, IP) 

C  PLOTXY 

A  PROGRAM  TO  PLOT  A  SIMPLE  REAL  VARIABLE  IN  A  SINGLE  PLANE 
INTENDED  AS  TOOL  FOR  ENGINEERING  PROGRAM  DEVELOPMENT.  IT  REQUIRES  A 
MINIMUM  KNOWLEDGE  OF  THE  GRAPHICS  SYSTEM  AT  THE  EXPENSE  OF  SOME 
'NICE  TO  HAVE  FEATURES,' NOTE A SLY  LABELING,  BUT  IS  EASY  TO  USE 
WITHOUT  HAVING  TO  LEARN  GRAPHICS. 

ThL  PROGRAM  CARO  MUST  BE  MOOIFIEO  TO  INCLUDE  GRAPHIC  FILE 
D£FINITION  AND  THE  FOLLOWING  ADDITIONAL  CAROSI 
PROGRAM  YOURS ( TAP  Ell, TAPE  12=97 /8?:> 

" 0 MMON  /CKT8UF/IPP,M8UFFt 1024) /CKXBUF /MDUMI 83) 

DATA  IPP/1S24/ 


THE  FORM  OF  THE  CALL  IS* 

CALL  PLOTXY(Y,NP,IR) 

WHERE  Y  IS  AN  ARRAY  OF  NP  OR  01 NA  TE  VALUES  TO  BE  P.OTTEO. 

THE  ARJ A  Y  WILL  BE  AUTOMATICALLY  SCAlEO  TO  FIT  IN  A  DE-AULT  SPACE 
PROVIDED  ON  THE  SCREEN.  WHEN  THE  PLOT  IS  COMPLETE,  THE  CJRSOR  WILL 
APPEAR,  PROVIDING  A  PAUSE.  ENTERING  ANY  CHARACTER  WILL  RETURN 
EXECUTION  n  THE  CALLING  PROGRAM  WITH' THE  CHARACTER  ERTEREU 
RETURNED  IN  IR. 

THAT  IS  ALL  Rt-QUIR EG  TO  OBTAIN  A  PLOT.  A  FEW  ADDITIONAL  FEATURES 
HAVE  BEEN  ADDED, BUT  MAY  BE  IGNORED.  THESE  FEATURES  ARE  ACCESSED  WITH 
THE  ROLL  OWING  CALLS) 

CALL  PLOTXY (Y,NP,IR,IP) 

WHERE  All  PARAMETERS  ARE  AS  BEFORE  AND  IP  SPECIFIES) 


3t  T 


2 


IP  = 


IDENTICAL  T0  PREVIOUS  CALL 

PLOT  Y,  BUT  USE  THE  SCALING  CALCULATED  BY  THE  PREVIOUS  CALL. 
SET  THE  PHYSICAL  AREA  ON  THE  SCREEN  TO  THE  ARIA  OUTLINED  BY 
TH'  NEXT  TWO  CURSOR  INPUTS. 

SET  THE  PHYSICAL  AREA  TO  THE  DtFAULT  AREA  PROVIDED 
SCALE  THr  OATA  IN  Y,  BUT  DO  NOT  PLOT. 

2  SET  PLOT  SIZE  ON  SCREEN 

0  AUTO  SCALE  AND  PLOT 

1  PLOT  TU  OL  C  SCALE 

3  SET  DEFAULT  SIZE 

4  AUTO  SCALE,  NO  PLOT 


C 


54 


APPENDIX  C 


SAMPLE  RUNS 
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4, 


•  ::  .  Ar. 


SAMPLE  COMPUTER  RJNf 

fu»ar  alpha  Input  in  louar  cast  f«r  clarity) 


Enter  J  F OR  INSTRUCTIONS,  DI  FOR  DEFAULT  UALUES 
N  Pax.  3?  UNUSED  CORE*  27722 

Enter  kevuord 


The  IPPEDANCE  op  tuo  thin  parallel,  strips  centered  uith  respect 

to  their  pidline  is  calculated.  . . 

enter  parapeter  or  function  kevuord  prop  belou  list  to  OPERATE. 

(USUALLY  THE  FIRST  TUO  CHAR  OP  DESCRIPTION) 

NO  KEVUORD  INDICATES  AN  OUTPUT  UARIABLE 

RESOLUTION  IS  OPTIPIZED  UHEN  IDCT*#,  PAV  IE  REDUCED  TO  IDCT 
DIOETS  IF  IDCT  IS  NON-ZERO.... 

DISPLAY  UILL  PRESENT  CURRENT  UALUES 


DESCRIPTION  KEVUORD  UNITS 
ELENENTS  N 

frequency  fp  hZ 

separation  SE  PETERS 

uiotmi  j:  peters 

widths  wl  peters 

ELE.uIDTHI  peters 

ELE.JIDTH2  PETERS 

thickness  n«  peters 

ski**  DPTHS  SK 

EFF. THICKS  PETERS 

RHO  RM  OMPSIPETER 

DIE.CC-nc  ED  RELATIVE 

POWER  FfcCT  AC  <SINi 

RESOLUTION  ::  DICETS 

CO  CO 

PLOT  PL 

PRINT  PP 

STOP  ST 

display 
instruct  I 

SAVE  TAPE2  SA 


ENTER  KEVUORD 
UARIABLE  UALUE 

elements  •  2.000 

FPE'jUiNC  y  •  :O0.0 


5Ea**FATICN' 

JIDT*U 

UIDTH2 

ELE.UIDTHI 

ELE.UIDTH2 

THICKNESS 

5KIN  DPTNS 


i.aoME-ea 

1.O000E-02 

1.0000E-02 

e. 

e. 

L.S**eE-®4 

;.O0e 


£FF,ThIC«:i' 

S"0 

r ie .cons 

P-JER  F4Ct 
rfSOLUTION 


1.72A1E-08 

1.000 

0. 


0. 


HZ 

PETERS 

PETERS 

PETERS 

PETERS 

PETERS 

PETERS 

PETERS 

OhPS* peter 

RELATIVE 

'SIN' 


ENTER  KEVUORD 


.02  SECONDS 


Mutual  solution  time* 

*IDGT  IN  INVERTING  PSI*  14 
UOLTAGE  DROP  in  BOTTOM  (UIDE)  TAPE* 

■  S  0. 

UOLTAGE  SOLUTION  TINE*  .01  SECONDS 

«IDGT  IN  P  SOLUTION*  13 
«IDGT  IN  I  SOLUTION*  14 
CURRENT  SOLUTION  TINE*  .01  SECONDS 

(DIFFERENCE  IN  CURRENTS*  0.  PERCENT 
THIS  RUN  REQUIRED  .04  SECONDS 


EFFECTIUE  INDUCTANCE*  G.2I32E-07  hENCRVS/N 
EFFECTIUE  RESISTANCE*  1.3793E-08  OHNS/N 
CHARACTERISTIC  INPEDANCE*  111.4  OHNS 
EFFECTIUE  CAPACITANCE*  1.770IE-11  FARAD0/M 
SHUNT  CONDUCTANCE*  0.  PHOS/P 

COMPLEX  IMPEDANCE*  701.0  OHMS  J  -770 
PROPAGATION  CONSTANT*  I.03S3E-O0  NEPERS/M  I 

7.S000E-0S  DI/M  S 


ENTER  KEYWORD 


LOWER  TAPE 

REAL  -J  ' I MAC) 

1 .S112E*01  4.7057E-01 

1.8109E40I  S.0015E-01 


CURRENTS 

UPPER  TAPE 
MAGNITUDE  REAL 
i.ani£«0i  *i.ni2E*0i  - 
l.tllSE»01  -1.01O9E+O1  -! 


ENTER  KEYWORD 
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1  OHMS 
NC0E-0C  RADIANS/M 
0O13E-04  DEGREES/M 


J  ' IMAG  >  MAGNITUDE 
I.70S7E-01  1.8110E-01 
I.001SE-01  1.8lllE«01 


|NTE*  N 

N  max.  37  UNUSED  CORE*  26098 
ENTER  kEvuORD 

enter  ualue 
9 

enter  keyword 

enter  10TN  STRIP  WIDTHS 
-••I  .MI 

enter  keyword 

EN?ER  VALUE 

.•MI 

enter  keyword 

Mutual  sour  ion  tine.  .26  seconds 

EFFECTIVE  TARE  THICKNESS  ADJUSTED 
tXDCT  IN  INVERTING  RSI*  l# 

VOLTAGE  DROP  IN  IOTTOH  (WIDE)  TAPE* 

■  SMMMMMIl  t.5«M4M47?92E-14 
VOLTAGE  SOLUTION  TIRE*  .34  SECONDS 


«i2gt  in  a  solution*  It 

•IDGT  IN  •  SOLUTION*  10 
CURRENT  SOLUTION  T«£. 


.•8  SECONDS 


•DIFFERENCE  IN  CURRENTS*  S.ISS»M434II2C-I2 
PERCENT 

THIS  RUN  REQUIRED  .61  SECONDS 

YOU  HAVE  2.214  CPU  SECONDS  LEFT 

AtOUT  ENOUGH  TINE  FOR  3  N OK  RUNS  LIKE  T»2$ 

EFFECTIVE  INDUCTANCE.  I.WS6E-07  NEnERyS'N 
EFFECTIVE  RESISTANCE.  17.16  OMMS/fl 
CHARACTERISTIC  IMPEDANCE*  32.SS  OhMS 
EFFECT lUE  CAPACITANCE*  I.0249E-10  FARADS'N 
SHUNT  CONDUCTANCE*  •.  MHOS/M 


COMPLEX  IMPEDANCE* 
PROPAGATION  CONSTANT* 


32.55  OHMS  J  -.4094 
.2636  NEPERS  'Pi  20.96 

2.29#  Df/M  1201. 


OHMS 

PAD IANS 'M 
DECREES /M 


ENTER  keyword 
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print 


CURRENTS 


LOWER  TAPE 

REAL 

-j  ciuag) 

HAGN1TUDE 

6.2920E-08 

1.2716E-04 

1.2731E-04 

2.6SSSE-07 

7.6123E-05 

7.6123E-05 

2.0448E-06 

7.8649E-0S 

7.867SE-0S 

1 •6197E-06 

7.6337E-05 

7.6354E-0S 

1.6723E-M 

7.5582E-05 

7.5601E-05 

1.6437E-08 

7.S040E-05 

7.S058E-0S 

l.SXK-N 

7.4718E-85 

7.4730E-05 

1.6295E-M 

7.4517E-05 

7.4S34E-0S 

1.S2C7E-M 

7.4426E-0* 

7.4444E-0S 

enter  keyword 

UPPER  tape 

REAL  -J  :inac:  "lAGNr'uDE 
-6.2920E-06  -i.2Tl6E-04  1.273;£-04 
-2.6S5SE-07  -7.6123E-0S  7.6;23£-es 
-2.e448E-06  -T.864SE-0S  7.867SE-05 
-1 . 6197E-06  -7.6337E-0S  7.63S4E-0S 

*1 .6723E-06  -7.5582E-05  7.S601E-0S 
-1.6437E-06  -7.S040E-0S  7.50S8E-05 
-1.6360E-W  -T.4718E-05  7.4738E-0S 
-1.629SE-06  -7.4S17E-05  7.4534E-0S 

-1 .6267E-06  -7 . 4426E-05  7.4444E-0S 
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plot 

PROGRAM  pauses  uith  cursor  after  each  plot,  enter  A  CHAR  tO  CONTINUE 
after  the  second  PLOT.  ENTRV  OF  THE  FIRST  LETTER  of  The  BELOU  DESCRIPTION 
WILL  CAUSE  THAT  RESULT i 
REPLOT  BOTH  TAPES 
FULL  WIDTH  SMALLER  TAPE  i  REPLOT 

PLOT  SIZE  CHANCE  (TWO  CURSOR  INPUTS  DEFINE  FRAME  SIZE) 


ENTER  KEYWORD 
•top 
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